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1 . INTRODUCTION 


Three scatterometer systems at frequencies of 0.4, 1.6, and 13.3 
GHz have been installed on a National Aeronautics and Space 
Administration (NASA) aircraft as additional microwave sensors 
for the Joint Soil Moisture Experiment (ref. 1), which is now the 
Agricultural Soil Moisture Estimation Program (ASME) . These 
scatterometers transmit radio frequency (RF) waves downward and 
measure the power of the waves backscattered from an agricultural 
scene. The final results of the operation are expressed in terms 
of the backscattering coefficient a* as a function of the inci- 
dent angle 0 of the RF waves. Because a® (0) at a given frequency, 
depends only on the scene parameters such as soil moisture con- 
tent, surface roughness, vegetative cover, and soil type, these 
scene parameters can be estimated (remotely sensed) by studying 
the magnitude and variations of a® with 0. 

The 13.3-GHz scatterometer system has been studied rather exten- 
sively in the past (refs. 2, 3, 4) . The results of recent study, 
which was based on the new laboratory measurements of the system 
parameters, was reported by Rosenkranz (ref. 5) . Using the same 
approach, an up-to-date analysis of the 1.6-GHz scatterometer 
system was reported by Wang (ref. 6). Reports on the 0.4-GHz 
system are limited to a general system description by Reid 
(ref. 7) and by the vendors (refs. 8, 9) . Therefore, because of 
component replacement and antenna redesign a current analysis on 
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the 0.4-GHz system is required in order to assess the performan- 
ces and limitations of the system. It is the objective of this 
report to provide an^ up-to-date ahaly sis of the 0.4-Gflz scatter- 
ometer system, , based on the most recent measurements of the 
system parameters. 


The. specif iCatiCTtr ponfigurati^ functional description of 

the ^Stem, are discussed ^in 2. This section also includes 

ihe signal frequency , spectrum processing" procedures , and 

the rece^ level at aircraft altitudes. The 

^sdatterom^^ syst€^ is giveii in section 3. In this 

section, the ba<%scatterihq^^^ c^ aigorithitns used to cal- - 

quiate a*? -are presehted. A- return signal power , which uses a V 
continuous sy Stein calib^^ apj^pach is traced; through th^ / 

receiver -channeis.-^ Resolution cell length requir^ents are also^V 
discussed rin section Section 4 inelud.es the antenna patterns . : 
datay-which are studied for the, configuration of main lobe, the ~ 
effects of side lobe, and the cross-pblarization levels. Section 
5 gives the. numerical results , including the es timateS: Of q® and^^ 
the effect Of cross-pdiarization levels. Section 6 qreseht's the 
major cohclusipris-Of this study and several suggestiPhs pn futhre: 
endeavors for studying the 4 -GHz scatterpmeter * ^ , 


2 . SYSTEM DESCRIPTION 


2.1 SYSTEM SPECIFICATION 

The 0.4-GHz scatterometer , which is flown on a Lyndon B. Johnson 
Space Center (JSC) aircraft, is an interrupted continuous -wave (ICW) 
Doppler radar system (refs. 7, 8, 9), designed to measure the back- 
scattering coefficient a® as a function of the angle of incidence 0 
for various types of terrain. The scatterometer transmits an 
ICW of 400.85 MHz radio frequency (RF) waves alternately from the 
horizontal and vertical antennas. The backscattered signals are 
received by both antennas and processed by individual receivers 
for each transmitted polarization. As a result of the forward 
motion of the aircraft, different Doppler frequency shifts are 
introduced in the return signals from different angles of inci- 
dence. The power of the return signal from a ground- resolution 
cell which is associated with a given incidence angle is finally 
retrieved by bandpass filtering at the corresponding Doppler 
frequency during data processing. 

The system specifications and" capabilities are given in table 
2-1. There are three manual models of operation depending on the 
altitude of the aircraft. There is also an automode operation 
in which the pulse width (PW) and pUlse repetition frequency 
(PRF) are automatically controlled by an altimeter. For ASME 
"application, only manual mode 1-5 is used. The 1-5 indicates 
aircraft altitudes from 1000 to 5000 feet (304.8 to 1524 meters). 


TABLE 2-1. - THE 0.4-GHz SCATTEROMETER SYSTEM CHARACTERISTICS 


Terras 

Transmit frequency 
Transmitter power output 

Modulation 
Pulse width (PW) 

Pulse repetition frequency 
(PRF) 

Transmitter/ receiver (T/R) 
isolation 

Number of antennas 
Isolation between antennas 
Antenna gain (two-way) 

Antenna beam width (two-way) 
Antenna alongtrack coverage 
Position of peak antenna gain 

Maximum skew angle 
Polarization 
Dynamic range 
Receiver sensitivity 
Receiver noise figure 


400.85 ± 1 MHz 

40.8 dBm (peak); 25.4 dBm 
(average) 

Interrupted continuous wave (ICW) 

Manual 1-5# 1.35 ± 0.1 us 
Manual 4-20, 6.9 ± 0.2 ys 
Manual 16-40, 29.5 ± 0.74 ys 

Manual 1-5, 42 ± 0.8 KHz 
Manual 4-20, 10.5 ± 0.25 KHz 
Manual 16-40, 4.9 ± 0.125 KHz 

120 dB during transmit 
Two, horizontal and vertical 
25 dB minimum 
~ 25 dB 
~ 18® 

0 to 60® cosecant squared beam 

44® (angle of incidence) horizon- 
tal antenna 

68® (angle of incidence) vertical 
antenna 

0.25® 

HH, HV, VH, W 
120 dB 

-148 dBm (6 Hz bandwidth) 

■V 20 dB 
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TABLE 2-1. — Concluded. 


Terms 

Signal-to-noise ratio 
(1500 ft. altitude) 

Calibration 

Receiver output 


> 10 dB 

Simultaneous with data processing 

Four fixed-receiver gate channels; 
one calibration channel; signal 
level - 0 to 5 volts 
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The scatterometer is a multiple data channel system. Multiple 
receivers and switches automatically separate the backscattered 
signal into the following data channels. 

a. Horizontally polarized return when transmitting a horizon- 
tally polarized wave using a fixed-receiver gate (HHFG) 

b. Vertically polarized return when transmitting a vertically 
polarized wave using a fixed-receiver gate (WPG) 

c. Vertically polarized return when transmitting a horizontally 
polarized wave using a fixed-receiver gate (HVFG) 

d. Horizontally polarized return when transmitting a vertically 
polarized wave using a fixed-receiver gate (VHPG) 

The scatterometer continuously calibrates the backscattered signal 
by injecting a calibrate signal into each channel and monitoring 
the transmitter power using a calibrate channel. The ratio of 
power return to power transmitted is obtained by measuring the 
ratio of power return to calibrate signal in the channel of 
interest and by measuring the ratio of transmitted power to cali- 
brate power in the calibrate channel. The ratio of a calibrate 
signal in the data channel to calibrate signal in the calibrate 
channel is a predetermined system constant. This calibration 
method eliminates uncertainties which/ otherwise, would be 
introduced because of transmitter output and receiver gain varia- 
tion. 
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2.2 SYSTEM CONFIGURATION^ 


The 0.4-GHz scatterometer was installed in the NASA C-130 air- 
craft. The horizontal and vertical antennas were mounted near 
the aircraft tail such that the plane of the antennas made a 22® 
angle with the horizontal plane. Figures 2-1 and 2-2 show the 
side and back views of the antennas, respectively. Each antenna 
consists of 24 dipoles. Every crosstrack array contains eight 
equally spaced dipoles; there are three crosstrack arrays equally 
spaced in the alongtrack plane to form the complete antenna system. 
The dipole orientation of the horizontal and vertical antenna 
arrays are orthogonal to each other so that either the vertically 
or the horizontally polarized, cosecant-squared beam can be 
generated. The antenna pattern was shaped to cover an illumina- 
tion area bounded by a crosstrack angular extent of approximately 
18® and an alongtrack angular extent ranging from nadir to approxi- 
mately 60® in the aft direction. 

The transmitter/receiver (T/R) unit and the tape recorder were 
placed in a rack inside the aircraft. Coaxial cables were used 
to connect the T/R unit input and output to the respective antenna 
port. 


2.3 PRINCIPLE OF OPERATIO N 

The operational principle of the 0.4-GHz scatterometer has been 
described in several reports (refs. 7,8,9). A brief description 
is presented here to facilitate understanding the signal analysis, 
which is discussed in section 3.2. 
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Figure 2-1. — Side view of C-130 showing coordinates for 
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Figure 2-3 shows the block diagram of the scatterometer system. 

(It is a slight modification of figure 2-1 in reference 7.) 

The functions of major components of the system are described 
in the following paragraphs. 

The transmitted signal is an interrupted continuous-wave (ICW) 
whose width and duty cycle are controlled by the synchronizer. 

The frequency of the transmitted signal is determined by a 
crystal controlled oscillator in the transmitter/modulator. In 
addition to providing the 20-watt peak 400.85 MHz ICW output sig- 
nal, the transmitter/modulator provides a 400.85 MHz signal to 
the first local oscillator (LO) source and a transmitter sample 
signal (approximately 50 dB down) to the calibrate receiver. 

The timing signal is generated by the synchronizer. These 
include the modulator trigger supplied to the transmitter /modu- 
lator and the switching signals to the RF and intermediate fre- 
quency (IF) switches. The pulse width (PW) and the pulse repe- 
tition frequency (PRF) are determined by the mode selected by 
the operator. In the manual mode, the PW and PRF are fixed by 
the selected altitude range (table 2-1) . 

The scatterometer timing waveforms are shown in figure 2-4 . The 
transmitter/modulator is turned on by the modulator trigger pulse 
and has a fixed 0.2ysec delay. The RF energy from the transmitter/ 
modulator is alternately switched between the horizontally polarized 
antenna and the vertically polarized antenna by the antenna switch 
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Figure 2-3.— The 0.4-GHz scatterometer block diagram 
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Figure 2-4.— Scatterometer timing waveforms. 
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toggle waveform. Both receiving channels are coupled to their 
respective antenna by the antenna switch waveform. An inherent 
delay is built into this waveform to ensure that the transmitter 
is turned off before the antennas are switched to their receivers. 
The first LO switching waveform controls the LO input to the 
first receiver mixers. The first LO excitation is removed during 
transmission to provide additional transmitter/receiver isolation. 
Further isolation is accomplished in the IP switches, where 
receiver gating is done. The signals received from the antennas 
are gated by the antenna switch and applied to a directional 
coupler which provides a means of injecting the calibrate signal. 
The gated signal is attenuated very little by the coupler (i.e., 
introduced by the insertion loss of the coupler) and the cali- 
brate signal is attenuated by 20 dB (i.e., the coupling value). 

There are three amplification stages in the receiver units; namely, 
the RF amplification, IF amplification, and the post IF amplifi- 
cation. The combined backscattered and calibrate signal from the 
directional coupler are first amplified by a 27-dB gain RF ampli- 
fier which has a noise figure of less than 4 dB. The amplified 
signals are then fed to a mixer and down converted to the 60 MHz 
IF amplification stage. The input signal to the IF amplifier 
are subjected to an adjustable gain control in 10 dB steps over 
a 50-dB range. The output signals from the 60-MHz IF amplifier 
is fed to the IF switch and multicoupler assembly, where they are 
gated and distributed to the proper post-IP amplifier. The 
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post-IF amplifier provides both manual and automatic gain control, 
each having a range of 40 dB. The 60-MHz signals are first con- 
verted down to 5 MHz using the 55-MHz second LO and then converted 
down again to 500 Hz using the 4.995-MHz third LO. The 500-Hz 
signals are fed to a low-frequency amplifier where a direct current 
(dc) offset is added to form an output compatible for tape record- 
ing. 

There are four data channels, one calibrate channel, and three 
telemetry signals in the scatterometer system output. Each of 
the four data channels occupies a bandwidth of 1 kHz and is main- 
tained at a nominal level of 0.9 V root mean square by an auto- 
matic gain control (AGO circuit. In addition, the alternating 
current (ac) output signal is superimposed on a dc offset voltage 
of 2.5 Vdc. The calibrate receiver channel also contains IF 
and post-IF amplifiers, and the output consists of two discrete 
data signals and three discrete signals for telemetry of control 
positions. The data signals are located at 500 Hz and 1 kHz for 
the transmitter sample and the calibrate signal, respectively. 

The telemetry signals are obtained from three voltage-controlled 
subcarrier oscillators designed for operation on bands 1, 2, and 
3 of the Inter-Range Instrumentation Group (IRIG) standard. The 
mode selector switch is assigned to band 1; the horizontal gain 
control, band 2; and the vertical gain control, band 3. 
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2.4 SIGNAL FREQUENCY SPECTRUM 


To investigate the waveform of transmitted signal and its power 
spectrums distribution, a signal frequency spectrum is discussed 
in this subsection. In figure 2-4, the RF output of 400.85 MHz 
CW can be represented by A cos 2irf t [fig. 2-5 (a)]. The modula- 

O 00 

tor trigger provides a unit pulse train, S(t-nT), with 

n=-“ 

period T and pulse width t as shown in figure 2-5 (b). The trans- 
mitter output is a combination of the two sources [fig. 2-5 (c)]. 
This signal output from the transmitter can be expressed by 

00 

f (t) = A COS ZiTf^t ^ S^t-nT) (2-1) 

N=-co 

where 

S(t-nT) = 1 for nT-^<t<nT+j 

= 0 elsewhere n= 0, ±1, ±2, ... 

A = amplitude of the ICW at the transmitter output, and 
f^ = frequency transmitted =• 400.85 MHz 


The frequency spectrum is obtained by applying Fourier tr^msform 
on f(t). It is shown in reference 10 that if f(t) = f^(t) f 2 (t), 
then F[f(t)J = F[fj^(t)l ® F[f 2 (t)], where F[ ] denotes the Fourier 
transform of the function in the bracket and 0 denotes the convolu 
tion integral of the two functions in the expression. Now let 


and 


f ^ (t) = A cos 2TTf^t 


(2-2) 
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( 2 - 2 ) 


f 2 (t) = 2 S(t-nT) 


n=-oo 


Then, the transform of fj^Ct) is 


F(f^ it)] = A / cos 27Tf te 32 irft 


I 

■ii: 


-j2u(f ^ f^)t ^ ^-j2,r(f + f^)t 


o' dt 


= ^ 6(f - f ) + 6{f + f ) 


(2-3) 


where <5(f - f^) and 6(f + f^) are the delta functions which 
represent unit impulse at f = f^ and f = -f^, respectively. 


The transform of f 2 (t) is 

-oo 00 

F[f2(t)] = I " nT)e"^^^^^ dt 

Loo 


(2-4) 


Observe that for n = 0 


f ”s(t)e^2’'“ dt = at = 

J-T/2 ^ 


and for n = K 


/ 

•/— C 


(2-5) 


S(t - KT)e 


-j2TTft 


-' — 00 
A 

^ -]2TTfKT / ^-j27Tff df 


/2 


sin TrfT _-j2irfKT 
Trf ® 


(2-6) 


2-15 


Thus 






(2-7) 


Therefore, the Fourier transform of f(t) is 


C 

— c 


F[f (t)] = J / [5(f - f^ - V) 


+ 6(f + f - V) sin_Z^ 

O ITT 


-j2TrnvT 


A T^(f “ fo)T ^ ^-j2Trn(f - f^)T 

2 7T(f - f ) Z-» 

o n=-“> 


A if(f + fQ>T ^ ^-j2-rrn(f + f„)T 

■^2 TT(f + f ) 


( 2 - 8 ) 


The frequency spectra, as derived in equations (2-1) through (2-8) , 

oo 

are shown sequentially in figure 2-6. The summation, ^ g-j2trnfT^ 

n=-oo 

gives unit impulse train spacing at PRF, f ^ [fig. 2-6 (a)]. 

si.n TrfT 

The shape of the function — is given in figure 2-6 (b) . 

The spectrum of F[f 2 (t)], [fig. 2-6 (c) ] , results from the multi- 
plication of figures 2-6 (a) and 2-6 (b). Finally, the convolution 
operation of F[f 2 (t)] with respect to F^fj^ (t) ] transfers the total 
spectrum to the center positions at f = f^ and at f = ”fQ. Figure 
2-6 (d) shows the spectrum of F[f(t }‘l Only the spectrum positioned 
at f = f^ will be of interest in the present discussion. 
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The spectrum of a transmitted signal with a 42-KHz PRF, a 1.35- 
Usec PW, and positioned at the 400-MHz carrier frequency is shown 
in figure 2-7. The spectrum of the received signal after a third 
down conversion is also shown. 


2.5 RECEIVER SENSITIVITY 

One way to estimate the sensitivity of the scatterometer system i 
by examining the signal-to-noise ratio at the receiver output. 

The power P^. reflected from a target and collected by a receiving 
antenna can be expressed in terms of the radar equation (ref. 11) 


as 


^r = 


(4tt) ^R^ 


(2-9) 


where 

P^ = power transmitted by the radar in watts 
= gain of the transmitting antenna in watts 

= gain of the receiving antenna in the direction of the target 
a = radar cross section of the target 
A = radar wavelengch 

R = distance between the radar and the target 

To compute the signal-to-noise ratio at the receiver output, the 
radar equation needs to include the effects of duty cycle and 
several loss factors; it is derived as follows. The antenna 
pattern factor can be expressed by; 
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where is the two-way antenna pattern factor. 


The radar cross section of the target can be expressed by; 

®(0)dA 




where a® (0) is the backscattering coefficient per unit area; the 
integration is taken over the illuminated region. 


The differential area dA is 

dA = R^d0d(J) (2-12) 

Substituting equations (2-10) through (2-12) into equation (2-9) 
yield 



The Doppler frequency at the angle of incidence (0^) and the air- 
craft ground velocity (V) is expressed as; 

«o 

and 



2V 

T 


cos 6 A6 


(2-15) 
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Since a wide aft beam antenna is being considered, fj^(G,(|)) is a 
slowly varying function of 0 and may be written as 



f^(0,(|>)d(() = B^F (0) 


(2-16) 


where B. is the two-way crosstrack effective antenna beam width, 
s' 


The power density in watts per hertz r received at the 

receiving antenna, is obtained by substituting equation (2-16) 
and the relation h = R cos 0 for equation (2-13) and by dividing 


by Afjj. 


W (0 ) = 
r o' 


P (0 ) 
r o 

Afr. 


P^G^F (e^)^^g°(6^)B^ cos 
2(4iT)\h^ 


(2-17) 


where h is the aircraft altitudes. 


Figures 2-8 and 2-9 show the resolution cell as defined by the 
Doppler frequency bandwidth and antenna beam width and the geom- 
etry of a scatterometer ground- resolution cell, respectively. 

The average power transmitted is P^^g = 

duty cycle of the transmitter. The duty cycle is expressed 
by : 

Di = T . £r 

where t and f^ are PW and PRF, respectively. 


(2-18) 

Of this transmitted 
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Figure 2 - 9 . ~ Geometry of a scatterometer ground resolution cell. 
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power, only D,p amount is usable; the remainder is redistributed 

as harmonics of the PRF. Therefore, the relative power contained 

2 

in the fundamental waveform (fig. 2-7) is 


The received signal power density becomes 


W (6) = 
rs 


cos e 


2(4TT)\h^L 


(2-19) 


m 


where is the miscellaneous loss, including the two-way antenna 

2 2 

coaxial cable loss and switching loss; G (0) = G^F(0) and 0^ have 
been generalized to 0. 


The noise power per hertz is expressed by: 

W = RTF D„ (2-20) 

R 

where F is the receiver noise figure, and is the receiver gated 
on duty cycle. 


The signal-to-noise power density ratio is 

W.^s PtG^(e)X^c;°(9)B^ cos 9 

'^n 2(4TT)^Vh^L KTFD„ 

m « 


( 2 - 21 ) 


The numerical values of the system parameters in equation (2-21) 
were evaluated using the measured data and were included in table 
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2-2 for the estimate of the signal-to-noise ratio. The aircraft 

altitudes and the speed were assumed to be approximately 457 m 

(1500 feet) and 77 m/sec (253 ft/sec) , respectively. The incidence 

angle was taken to be 60®/ where low signal return was expected. 

The a® value at a 60®-incidence angle was taken from figure 3-3; 

the calm sea value was used for low signal return. A duty cycle 

of 0.0567 for manual mode 1-5 was computed from equation (2-18). 

All of the remaining parameters appeared in equation (2-20) are 

self-explanatory. The gain of the receiver amplifier was not 

entered in the computation because it applied equally to signal 

and noise and was eliminated when the signal-to-noise ratio was 

taken. In calculating the receiver noise level, the ambient 

temperature (T) was assumed to be approximately 290 K. The 

“23 

Boltzmann constant (K) is 1.38 x 10 J/K. The 5.4-dB noise 
was obtained from the report by the Emerson Electric Company 
(ref. 9) . The signal power density received is the sum (in dB) 
of all items given in equation (2-19) and is -131.9 dBm/Hz. The 
receiver noise power density is -173.1 dBm/Hz, which results in 
a signal-to-noise ratio of approximately 41 dB. However, the 
signal-to-noise ratio of the tape recorder system as reported iii 
reference 6 is approximately 40 dB, which represents the best 
value attainable with the 0.4-GHz scatterometer and data recording 
system. ^ 
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TABLE 2-2. - THE 0.4-GHz SCATTEROMETER SIGNAL-TO-NOISE RATIO 


Parameters 

Power 

Level 

Transmitter output power P^ 

40.8 

dBm 

G^(0) (W gain at 0 = 60«) 

19.3 

dB 

B, (at 0 = 60®, W beam width is 23.3®) 
<P 

-3.9 

dB 

(wavelength - 0.75 m) 

-3.7 

dB 

o®(0) (0=60®) 

-40.0 

dB 

cos 0 (0 = 60®) 

-3.0 

dB 

l/[2(4Tr)^] 

-36.0 

dB 

1/V (aircraft velocity - 77 m/sec) 

-18.9 

dB 

1/h^ (aircraft altitudes - 457 m) 

-53.2 

dB 

1/L (coaxial cable loss) 

TC[ 

-8.4 

dB 

2 

D,p (duty cycle - 0.0567) 

-24.9 

dB 

Received signal power density 

-131.9 

dBm/Hz 

KT 

-174 

dBm/Hz 

Receiver noise figure (F) 

5.4 

dB 

Receiver gated on duty cycle (Dj^) 

-4.5 

dB 

Receiver noise power density 

-173.1 

dBm/Hz 

Signal-to-noise ratio 

41.2 

dB 
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2.6 DATA PROCESSING 


The output of the scatterometer is recorded by a Mincom Series 
110 Recorder/Reproducer System, Wideband Group I.^ The specifi- 
cations of the recorder for the three scatterometers can be found 
in reference 6. The channel number assignments for various data 
output of the scatterometers are shown in figure 2-10. The 0.4- 
GHz scatterometer outputs are recorded on five adjacent data 
channels (4, 6, 10, 12, and 14) of the same tape. 

The data-taking time is recorded through the IRIG Standard Time 
Code and the Aircraft Data Acquisitions System (ADAS) in channels 
9. and 2, respectively. In addition to the time information, ADAS 
acquires the aircraft parameters such as altitude, speed, roll, 
pitch, and drift angles. 

The time and aircraft parameters provided by ADAS are recorded 
sequentially in pulse code modulation (PCM) format and frequency 
modulated at 225 KHz. All of this information is derived during 
data reduction and is used in data processing for the derivation 
of a® as shown in figure 2-11. 

Figure 2-11 is a brief data processing block diagram of the Texas 
A&M University software processor which processes the three 
scatterometer analog data tapes. In addition to the recorded 

T " — 

This model is used for identification purposes only and does not 
constitute a recommendation. 
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Figure 2-10.— Channel assignment of signal outputs for three scatterometers 
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Figure 2-11.— Data processing block diagram 




signals and ADAS data in the tape, the system calibration con- 
stants, the resolution cell length, and the angle of incidence 
are required as inputs. Digital filtering is employed in the 
processor to obtain the signal spectrum at the required center 
frequency and bandwidth which corresponds to the specific inci- 
dence angle and area, respectively. At the same time, the compu 
tations of antenna gain and area, with pitch, roll and drift 
corrections, result in the computation of a° versus 0^. The 
output of o® is tabulated with respect to incidence angle and 
flight time. 
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3. SYSTEM ANALYSIS 


3.1 BACKSCATTERING COEFFICIENT ALGORITHMS 

The power density W^(0) at the receiving antenna is given in 
equation (2-17). Expressed in terms of a® (0) / equation (2-17) 
becomes 


2 (4tt) \h^W (0) 

a®(0) = — (3-1) 

cos 0 G'^(0)B.P. 


Since W^(0) = P^(0)/Afj^, equation (3-1) can also be written a:.: 


c®(0) 


2(47r)\h^P^(0) 

cos 0 G^(0)B. Af^P^ 
<P D t 


(3-2) 


where Af^^ is related to the ground cell length Ap^ in the follow- 
ing form (ref. 11) . 


, £ 2V cos 0 * ^ 
= Xh 


(3-3) 


Equation (3-2) is the simple, direct, and most used approach of 
determining 0 ® as a function of incidence angle, 0. As derived 
in the previous section, it associates a ground resolution cell, 
A(0), with a single average value of a® (0) (fig. 2-8) and then 
expresses o® (0) in terms of all measurable parameters. 
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The a® derived from this approach contains two sources of 
error, both resulting from the approximation of a rectangular 
area (fig. 2-9) for an actual Doppler cell defined by two iso- 
dopplers (fig. 2-8) . The first source of error comes from 
the assignment of a single value of a® for a given rectangular 
ground cell. Actually, a° is observed to vary with incidence 
angle 0, and the contribution to the observed signal return 
within a given bandwidth, Afj^, is the sum of the backscattered 
signals from all the small area elements with different 0^’s in 
the Doppler cell. The second source of error is the results of 
crudeness of the area approximation itself. Both of the error 
sources depend on the scatterometer antenna beam width. For 

the 0.4-GHz scatterometer, the antenna beam width, B,, varies 

<P 

from approximately 14° to 25° over the range of 0, from 0° to 
60° (fig. 4-5) . In this case, it is desirable to explore the 
amount of error introduced in the a° derivation by using equa- 
tion (3-2) . 

For deriving received power, P^., without error-introduced approxi- 
mation, refer to the radar equation 


P = 
r 


P. G. G oX 
t t r 

(4tt) 


(3-4) 


The radar cross section of the target a is the integration of 
backscattering coefficient a° through the differential area 
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element, dA, and is expressed as: 

do = o®dA (3-5) 

In terms of equation (3-5) , the radar equation becomes 


f P. G. G 0 
t t r 

(47T) \ 


T 


dA 


(3-6) 


The power, Pj^(O), observed at a given central Doppler frequency, 
fjj(G), and within a bandwidth, Afp(0), is uniquely related to a 
ground cell. A, as shown in figure 2-8. The integration in 
equation (3-6) is computed by numerical summation because only 
a discrete value of measured antenna gain pattern is available 
and the ground cell. A, is further divided into N elements. 
Equation (3-6) is then expressed as: 


P_(6) = 


!iL 

(4tt) 


i=l 


,)G3^(0^)o^(0^) 
R^(0^) 


AA. 


(3-7) 


The power density, W (0) , is simply 


P„(6) 

«r<«> = 


4fp(e) (4n)^ ( 0 ^) ^ 


(3-8) 
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Examining equations (3-7) and (3-8) , it is obvious that no 
explicit expression for a® (0) is available for the exact deriva- 
tion which makes rqu<ition (3-7) or (3-8) very difficult to use 
for data reduction purposes. But it can be used for comparison 
to evaluate the error induced in equations (3-1) or (3-2) when 
large antenna beam width and Doppler bandwidth are used. 

To derive the backscattering coefficient, a®, as a function of 
incidence angle, 0, from equation (3-8) requires an iteration 
process. First, a functional dependence of a® on 0 has to be 
assumed. Then the power density, W^(0), of the backscattered 
signal is calculated from equation (3-8) for a series of 0 and 
A0. The calculated W^(0) is compared with the power spectrum 
which is derived from the observations. If the computed and 
observed power spectra do not match well over the incidence 
angles of interest, another form of a® (0) is assximed and the 
procedure is repeated. The final functional dependence of a® 
on 0 is attained when the best match between the computed and 
the observed power spectra is obtained. 

Equations (3-1), (3-2), (3-7), and (3-8) apply to each of the 

four polarization combinations as described in section 2. Thus, 

P / P._f G , and can be in either the horizontal or the ver- 
r t r t 

tical polarization state. Depending on the polarization state 
of (or G^) and P^ (or G^) , four different backscattering 
coefficients can be derived from the operation of the scatterometer ; 
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namely, ^The subscripts indicate 

the polarization state of P^, P^, G^, G^, and o®, which are not 
explicitly written in the expression above.) The power received, 
P^, is the power of the backscattered signal at the receiving 
antenna output port. As the signal passes through the receiver 
data channels, it experiences losses as well as gains because of 
the presence of many subsystems. These factors will have to be 
considered in order to correctly interpret the observed data. 

Both approaches for the determination of a® versus 0 relationship 
[namely, the a® values from equations (3-2) and (3-7)] are 
attempted in section 5, where the numerical calculation is per- 
formed. A comparison is made on the results from both approaches. 

The backscattered signal in the receiver data channel and its 
continuous calibration technique is analyzed further, and the 
backscattering coefficient, a®, is discussed in the following 
sections. The antenna gain pattern, side lobe level, and cross 
coupling effects are discussed in detail in section 4. 

3.2 RECEIVER SIGNAL ANALYSIS - CALIBRATION CONSTANTS APPROACH 
The desired data in equations (3-2) and (3-7) are the relation- 
ship between P^ and P^. If this relationship is established at 
a common point in the system, then combined with antenna pattern, 
aircraft velocity, and aircraft altitude data, it is possible 
to obtain the desired backscattering coefficient. The antenna 
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gain data must be referenced to the common point selected in 
the system. In the 0.4-GHz scatterometer system, calibrate 
signals are generated, which serve as a reference in the 
receiver channel . These calibrate signals provide information 
in the recorded data which can be used to eliminate uncertainties 
which would otherwise be introduced because of transmitter out- 
put and receiver gain variations. 

The various signals in the 0.4-GHz scatterometer system are 
indicated in figures 3-1 and 3-2. Figure 3-1 shows the rela- 
tionship between the received signals and the calibration sig- 
nals. As indicated, the desired relationship is the ratio of 
to P^. Referring to figure 3-2, the transmitted power, P^^ 

(at the horizontal antenna port) , is equal to the peak trans- 
mitter power Pq less the attenuation through the antenna switch 
and the RF filter. 


^TH *^THl^TH2^0 


(3-10) 


A sample of the transmitted power is fed to the calibrate 

receiver through a directional coupler. The coupling is desig- 
nated as K^ 2 * Thus, the transmitted power sample in the cali- 
brate receiver has a value of calibrate signal 

which is supplied to the calibrate receiver is designated C^; 
the value at the front end of the calibrate receiver is 
where is the constant introduced by the insertion loss of the 
directional coupler. 
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Figure 3-1.— Received signal calibration 












Figure 3-2.— Transmitted signal calibration 



The receivers used in the 0.4-GHz scatterometer pass only the 
central spectral line. Thus, the transmitter sample signal 
from the calibrate receiver is modified by the square of the 
transmitter duty cycle. > 

2 

Transmitter sample = G (3-10) 

C T2 T± TO 

where G is the gain of the calibrate receiver, and D. is the ~ 

transmitter duty cycle. The calibrate signal from the calibrate 

receiver is equal to G K C . 

^ c cc c 

The received signal is processed in four separate receivers. 
Figure 3-1 shows the relationships between the received signals 
and the calibration signals in the receiver channels. The 
relationship of the horizontal receiver will be discussed here. 
Similar relationships are valid for all of the fixed-gate 
receivers . 

The received signal, is attenuated by the RF filter, the 

Kn 

antenna switch, and the directional coupler. The received sig- 
nal in the horizontal receiver is given by ^rjii^PU2^RH3^RH* 

The calibrate and received signals from the horizontal receiver 
are: 

Calibrate signal » ®rh^CH1^CH3^CH4^R^R1 (3-11) 

Received signal » ®rh^rh1^RH2^RH3(“i) ^RH (3-12) 
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For equations (3-11) and (3-12 ) , it is assumed that received 

signal eclipsing did not occur because of IF or LO switching. 

The duty cycle of the received signal in the fixed-gate channels 
1 2 . 

is since a particular receiver channel is active on alternate 
' pulse intervals. 

-Figure 1-rl indicates that two calibrate signal levels are avail- 
^j^alS^e receiver. These signals are alternately 

switched; one iS used during the active time of the low-gain 
channel and the other is used during the active time of the 
high-gain channel. 


The following relationship can now be established. 

Transmitter sample _ * _ ^c^T2^Tl^T^o 
Calibrate signal G K C 

w w w ^ 


(3-13) 


The gain term cancels out and the relationship becomes 


^T1*^T2^T^o 

* = --r-e — 

cc c 


(3-14) 


For the horizontal receiver channel under consideration the 
relationship is: 


Received signal _ „ _ ^RH^RHl^RH2^RH3 T? ) ^RH 
Calibrate signal “ “ I ~ ~ ~ 17 _ 

^Rh'^CH1^CK3^CH4^RH^R1 


(3-15) 
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Removing the gain terms yield: 


_ ^RHl^RH2^RH3^T^RH 
^^CH1^CH3^CH4°RH^R1 


(3-16) 


The values of A and B are measured quantities available from 
the recorded data. 


The relationships given in equations (3-14 ) , (3-16) / and (3-9) 
can be combined to give; 


RH 

*TH 


^^ ^CH1^CH3^CH4^RhS i 
^RH1*^RH2^RH3 


B 




^cc^THl^TH2 
^Tl^T2 i 


iCcA 


Let 


/ ^^CH1^CH3^CH4 \ ^ 

\ ^RH1^RH2^RH3 / 

/^CC^TH1^TH2\ _ ^ 

\ ^T1^T2 ) ^ 2HH 


(3-17) 


(3-13) 


Then 


^RH ^ / ^IHH * ^3Hh \ q2 , B 

^TH \ ^2HH / R A (3-19) 
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Each receiver has a value for and K^. Dj^ is established 

for each receiver and is a function of mode and range of opera- 
tion. 

The calibration constants are determined by measuring the indi- 
vidual components to establish values for and K 2 . The value 
for requires a direct measurement of the power levels of 
C„, and which are fixed power levels for all operating 
conditions. The quantity D_, which is variable with the mode 
and raiige of operation, is established by measuring the calibrate 
signal level. In these measurements, the AGC is inoperative so 
that the receiver gain does not vary from measurement to measure- 
ment. 

The system calibration constants are expressed by: 

y _ /^1HH^3Hh\ 

\ K2JJJJ j R (3-20) 

When K^, K 2 / and are obtained by direct measurement, the 
values of the calibration constants (K) can be computed from 
equation (3-20) . The values are listed in table 3-1 for differ- 
ent polarizations, gain settings, and modes of operation. 

The calibration constants listed in table 3-1 are the most 
recent data which reflect the changes introduced by the replace- 
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ment of parts and the addition of attenuators in the calibrate 

source. The values originally given by the Emerson Electric 

Company (ref. 8) were corrected by -22.7 dB for linear channels 

and -24.5 dB for cross channels. Corrections were made because 

of the installations of spare transmitter coupler and the output 

power divider of the calibrate channel with 20 dB pads added. 

Linear +22.7 dB and then dropped to +3.8 dB. Cross 

Cj^^/C^ was +11.62 dB and then dropped to -9.1 dB. The two-way 

antenna coaxial cable loss (l ) of 8.4 dB was also measured in 

s 

Building 15 of JSC. 

The calibration channel coding is shown in table 3-2 . From the 
frequency spectrum given in the calibrate channel, the values 
of the calibrate signal and the transmitter sample can be ob- 
tained along with the receiver gain position and the mode of 
operation. The calibration constants can be selected from table 
3-1 when the receiver gain position and the mode of operation 
are known. 

3.3 MEASURED SYSTEM PARAMETERS 

After establishing the ratio of to by using the calibra- 
tion constant for a specific polarization state, the gain setting 
of the receiver, and the measured antenna cable loss L , it is 

5 

now possible to put all the parameters together and derive the 
final forms of the backscattering coefficient algorithms. The 
expression for equation (3-2) becomes 
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TABLE 3-1.- CALIBRATION CONSTANTS 
[Values in negative dB] 


Manual (1-5) mode, fixed gate 


Channel 

Gain position 

1 

2 

3 

4 

5 

6 

W 

112.5 

123.0 

132.8 

142.3 

152.3 

162.5 

HH 

112.8 

122.8 

132.8 

142.7 

152.1 

161.8 

VH 

125.7 

135.7 

145.7 

155.6 

165.0 

174.7 

HV 

125.4 

135.9 

145.7 

155.2 

165.2 

175.4 


TABLE 3-2.- CALIBRATION CHANNEL CODING 



Code frequency, Hz 


Manual (1-5) 
Manual (4-20) 


Horizontal receiver 
gain position 


Vertical receiver 
gain position 


Transmitter sample 
Calibrate signal 



500 

1000 
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(3-21) 


0(0) ^ 2(4tt)\A(B/A) 

cos 0 G^(0)B. Af_L„ 

9 D S 


where 

B = the ratio of received signal to calibrate signal in the 
receiving channel 

A = the ratio of transmitter sample to calibrate signal in the 
calibrate receiver 

L = two-way antenna coaxial cable loss (-8.4 dB) 
s 

K = calibration constant for specific receiving polarizations 
and gain settings 


ilCXVC 


The other parameters in the equation 
Similarly, the expression for equation (3-7) becomes 

(B/A) _ 


< 


Af, 


N G. (0.)G (0.)a°(0.) 
t 1 r 1 1 

K Afp(4n)* ^ R*(0^) 




AA. 


(3-22) 


Note that in equation (3-22) the ratio of B over A is the relative 
receiver output power with respect to the transmitter sample power 
in the Doppler frequency band (Af jj) , which corresponds to the 
signal return from the ground cell defined by two isodopplers 
(fjj + j Afjj and j Af^^) for a specific angle of incidence (0) 

which generates fj^. 
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The values of K and could change for different missions. The 

calibration constant K controls the power level of the calibrate 

signal, which is subject to the limitation of the dynamic range 

of the tape recorder. For a given surface of backscattering 

and aircraft altitude, K could be adjusted for an appropriate 

calibrate signal level using different values of attenuation 

pad to bring about the best signal-to-noise ratio for the data. 

The value of depends on the length of the cables between the 

transmitter/receiver unit and the antenna, which can be changed 

from mission to mission. The values of K and L , measured for 

s 

the functional check flight of Mission 347 will be used for the 
nxamerical calculations in section 5. 

3.4 RESOLUTION CELL LENGTH REQUIREMENT 

In the software processing algorithm, the resolution cell length 
L is assumed constant over all incidence angles and is expressed 
by: 


L = Ap^ + Ap^ (3-23) 

where Ap^ is the resolution length corresponding to the Doppler 
frequency bandwidth Afj^ as given in equation (3-3) , and Ap^ is 
the resolution length corresponding to the integration time. 

V X p 

Ap^ = — _s (3-24) 
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In equation (3-24) , fg is the sampling frequency of the digitizer, 
and Pg is the number of sample points for one cell as used in the 
program. For constant L and Ap^, Ap^ will be the same for all 
incidence angles; and from equation (3-3) , it is easily seen that 
Doppler bandwidth Afj^ will decrease rapidly as the incidence angle 
increased from 0 to a larger value, say 60®. 

Because scatterometers are, by definition, calibrated devices used 
in determining scattering information on an absolute scale, fading 
poses particular problems. Fading is caused by the constructive 
and destructive interference of the coherent signal compoii£;nts 
reflected by individual scatterers within a resolution cell. 

Making a single, instantaneous measurement of scattered power will 
introduce errors since there are no means by which we can determine 
how much fading is taking place. Thus, it is necessary to make 
a number of measurements which can be averaged to reduce fading 
effects. The signal must be observed long enough so that sufficient 
independent samples may be collected when averaging over a time 
period or spatial distance to reduce fading effects and thereby 
obtain meaningful values of a®. Independent samples may arise 
from two causes: (1) independence within the Doppler bandwidth 

while illuminating a ground path from different angles and (2) a 
combination of return signals measured from completely separate 
ground patches . The number of independent samples occurring 
during the time the scatterometer passes a resolution cell length 
of L and a Doppler bandwidth of Afjj is given (ref. 11) by; 
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2L Ap£ cos' 


G 


Xh 


(3-25) 


The relative standard deviation of the scattering coefficient is 
expressed by: 

a _ i _ / All 

^ l2L Ap^ cos^ 8 


^ 1/2 


( 2 - 2 ^^, 


Because of the long wavelength of 75 centimeters (29.5 inches) of 
the 0.4-GHz scatterometer , a tradeoff must be made between a 
large Doppler bandwidth, which improves the data accuracy, and a 
small ground cell, which improves the length resolution of the 
data. Table 3-3 gives the number of independent samples and the 
standard deviation values for the 75-meter (246-foot) cell length, 
the 2500-Hz sample frequency, 1024 sample points, 77-meter 
(252.5 feet) per second aircraft ground speed, and 442 meters 
(1450 feet) altitude. Also included in the table are the center 
frequency and the Doppler bandwidth at incidence angles ranging 
from 0® to 60®. 


There are less than 10 independent samples in the resolution cell 
when the incidence angle becomes larger than 40®. To attain more 
independent samples and thereby increase the accuracy, the cell 
length must be increased, especially at large incidence angles. 
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TABLE 3-3.- INDEPENDENT SAMPLE T^D STANDARD DEVIATION 
FOR CONSTANT CELL LENGTH OF 75 METERS (246 FEET) 


Incident 

angle, 

0 

Central 

Doppler 

frequency, 

^D 

Doppler 

frequency 

bandwidth. 

No . of 

individual 

sample, 

N. 

1 

Standard 

deviation, 

a/u 

5 

17.9 

20.0 

19.4 

0.227 

10 * 

35.7 

19.3 

18.7 

.231 

15 

53.2 

18.2 

17.7 

.238 

20 

70.3 

16.8 

16.3 

.248 

30 

102.8 

13.1 

12.7 

.280 

40 

132.2 

9.1 

8.8 

.330 

50 

157.5 

5.4 

5.2 

.438 

60 

178.1 

2.5 

2.5 

.638 
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3.5 THE FUNCTIONAL DEPENDENCE OF P-BAND 0° ON 0 


The ordinary approach of deriving the backscattering coefficient 
0 ® by means of a Doppler shift frequency in the received signal 
frequency spectrum, such as the 0.4-GHz scatterometer , is to use 
the backscattering coefficient algorithm given by equation (3-21) . 
The crosstrack dimension of a resolution cell is determined by 
the half-power beam width of the antenna, and the alongtrack 
dimension of a resolution cell is determined by the bandwidth 
Afp of a filter, centered at the Doppler frequency fj^, which 
corresponds to a selected angle of incidence. It was pointed out 
previously that as B^ and Af^ become larger, (especially when 
is too large to ignore for the 0.4-GHz scatterometer) , a signifi- 
cant amount of error may be introduced in the derivation of a®. 

To explore the amount of error introduced and to present an ade- 
quate standard for comparison, equation (3-22) is presented as the 
second form of the backscattering coefficient algorithm. 

The derivation of o® in equation (3-22) required an assumed 
functional dependence of o® on 0. To see how the shape of the 
a® versus 0 curve affects the error in a® estimation in equation 
(3-21) , the measured values of a® from both smooth and rough 
surfaces of backscattering are considered. For the smooth sur- 
face, the measurements by Guinard (ref. 12) at 0.428-GHz and 
for the HH polarization over a calm sea are adopted. For the 
rough surface, data presented by Dickey (ref. 13) at 0.4-GHz 
and for the HH polarization are used. 
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The results from both of these measurements are shown in figure 
3-3. Roughness data are not available for the terrain measure- 
ment as presented by Dickey (ref. 13), However, it is not impor- 
tant in the present case where the primary interest is to have 
the functional forms of a® dependence on 0 for the numerical cal- 
culations of equation (3-22) . 

A quadratic regression was performed on the data from the measure 
ments over the calm sea. The result is given by: 

a® = 5.71 X io"^0^ - 0.9710 - 2.85 (3-27) 

The land measurements data (ref. 13) covered dry and irrigated 
conditions, with averaging over the entire test site. The result 
of the dry field data curve, with second-order polynomial approxi 
mation and smoothing, is expressed by 

a® = 8.44 X 1O“^0^ - 1.010 + 9.85 (3-28) 

In both equations (3-27) and (3-28) , the unit of a® is in dB and 
0 degree. Both of these expressions are used in the numerical 
calculations in section 5.2, where the two forms of the back- 
scattering coefficient algorithms are compared. 
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Figure 3-3. The functional dependence of a® on 0. 
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4 . ANTENNA PATTERNS 


4.1 MAIN LOBE 

The antenna pattern data, both horizontal and vertical, were 
measured by the Physical Science Laboratory of New Mexico State 
University for every 2® interval in both the alongtrack and the 
crosstrack directions. The measurements were made with a set of 
two identical antennas without using a mockup to simulate the 
actual antenna system which is mounted on the C-130 aircraft. 
Consequently, the antenna pattern used in the present analysis 
might deviate from the real ones. The effects resulting from 
this deviation on the measurements of the backscattering coeffi- 
cient o® is not known. Calibration of the airborne system 
against the ground-based system may be necessary in order to 
determine the seriousness of this effect. 

Since each of the horizontal and vertical antenna is used for 
the transmission and reception of the RF powers, the same antenna 
pattern is generated for the transmitting and receiving modes of 
operation. The horizontally transmitted and vertically received 
antenna pattern (two-way HV polarization) will be the same as 
the vertically transmitted and horizontally received antenna 
pattern (two-way VH polarization) . Two-way gain of the HH and 
W polarizations is the one-way gain of the corresponding hori- 
zontal and vertical antenna multiplied by two in dB. To see the 
spatial variations of the antenna patterns, the contour maps of 
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the two-way gain of HH, W, and HV (or VH) polarizations are 
displaced in figures 4-1, 4-2, and 4-3, respectively. Each of 
the three maps covers the angular ranges of (f) from -40° to 40° 
and of 0 from 0° to 60°. The direction of nadir is located at 
0=0 and (p = 0. From these figures, it is clear that the RF 
powers of transmission and reception are mostly limited for 
in the region -40° to 40° (although the side-lobe level appears 
rather high at 0 == ±30°) and for 0 from 0° to 90°. The pattern 
for all polarizations (figs. 4-1, 4-2, and 4-3) displays an 
increase in gain with incidence angle and a rather high side lobe 
around nadir. 

The variations of the alongtraclc, two-way gains with the angles 
of incidence in the aft direction for all three polarizations are 
shown in figure 4-4. It is clear that the maximum of the two- 
way gains for HH and W polarizations appears at 0 = 66° and 
0 = 44°, respectively. On the other hand, the maximum of the 
two-way gains for HV (or VH) polarization appears at 0 =; 52° 
which is about half way between the maximum of HH and W polariza- 
tions. 

From equations (3-2) , it is observed that the bac)cscattered 

power, Pj., is proportional to G^(0), a°, and R Since both a° 

-4 

and R generally decrease with increasing 0, P^. would be small 

at large angles of incidence. For W polarization, the decrease 

2 

in the two-way antenna gain G (0) for 0 > 44° further depressed 
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Figure 4-1.— The two-way gain contour of HH polarization. 
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Figure 4-5.— The variation of the 3-dB becun width 
with angle of incidence. 
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is the results of the existence of the rather high side-lobe level 
in this region which generates the pattern broadening effect 
^figs. 4-6 and 4-7) . As expected for HV (or VH) polarization, 
the variation of the antenna beam width, with respect to the 
incidence angle, stays half way between HH and W polarizations. 

In general, the beam width variation for W polarization is from 
14® to 24®; for HH polarization, it is from 17® to 24®. From 
this, it can be stated that W polarization has better cross- 
track resolutions in the region of 0 from 12® to 35®. 


The second method of defining a beam width is given by: 



y*^^GQF(0) f (6,4»)d(j) 

GqF(0) 


(4-1) 


where and ^2 crosstrack angles where the null points 

of the main- lobe antenna gain occur. The function f(0,(()) defines 

2 

the crosstrack gain variation and G^F(0), which is expressed in 
numerical value and not in dB and has the same meaning as before. 
The results of the calculations based on equation (4-1) are 
tabulated in table 4-1 for incidence angles ranging from 0® to 
60® at 2® intervals. The corresponding 3-dB beam widths are also 
included in the table for comparison. It is clear from the 
table that for all angles of incidence and all polarizations the 
values of the beam width as derived from equation (4-1) are larger 
than the corresponding ones obtained from the 3-dB method. On 
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TABLE 4.1- COMPUTED TWO-WAY ANTENNA BEAM WIDTH IN DEGREE 


Incidence 
angle 0, 
degree 

HH 

Method 1 
W 

HV,VH 

: HI! 

Method 

vy 

2 

' ^ hv;vk 

0 

17.2 

20.7 

18.6 

17.8 

22.4 

19.3 

2 - 

16. 5 

20.7 

- 18.4 

17.7 

22.0 

19.6 

4 

17.3 

19.8 

18.5 , 


21.6 

19.3 

6 

16.5 

18.0 _ 

17.2 

n\7 

■ ,20.1 

17.9 

8 

17.2 

16.7 

16; 8 

18.0 

• 19.0 

: 17.8 

10 

18.0 

16.5 

17.2 

18.7 

18.3 

18.2 

12 

18 ;0 

14.7 . 

- 16.2 

19.1 

17.3 

17.5 

14 

18.0 

1;4>8 

16.3 

19.0 

16.3 

17.4 V 

16 

18.0 - 

14.8 

16.5 

18.9 

16.1 

17.2 - 

18 

18.5 

14.0 

16.0 

18.8 

15.9 

16.9 

20 

17.2 

14.5 

15.4 

18.1 

15.6 

16.3 

22 

17.7 

14.0 

15.6 

18.6 

15.3 

16.6 

24 

16 . 7 

14.0 

15.2 

18.4 

15.0 

16.5 

^ 26 

17.0 

14.5 

15.6 

17.9 

15.4 

16.3 

28 

17.3 

14.5 

15.7 

18.3 

16.2 

16.7 

30 

17.3 

16.0 

16.5 

18.3 

17.1 

18.4 

32 

16.7 

14.5 

15.3 

17.9 

15.8 

16.4 

34 

16.7 

14.7 

16.2 

18.0 

16.4 

17.1 

36 

16.3 

16.0 

16.2 

17.7 

17.6 

17.5 

38 

16.7 

15.3 

16.4 

17.4 

17.0 

17.5 

40 

17.0 

16.0 

17.0 

17.8 

16.6 

17.1 

42 

18.0 

16.0 

16.9 

19.7 

17.0 

18.0 

44 

18.0 

16.7 

17.3 

19.0 

17.6 

18.2 

46 

19.0 

17.0 

17.9 

19.7 

18.4 

18.8 

48 

19.0 

18.7 

18.5 

20.0 

19.6 

19.5 

50 

19.0 

19.3 

21.3 

20.2 

20.8 

20.4 

52 

20.3 

19.7 

21.7 

22.1 

21.0 

21.5 

54 

22.0 

19.3 

20.5 

23.1 

21.2 

22.0 

56 

22.5 

20.7 

22.0 

24.3 

22.6 

23.4 

58 

22.0 

22.0 

22.4 

23.2 

23.8 

23.2 

60 

24.0 

23.3 

24.2 

25.5 

25.0 

25.2 
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the average, the difference in beam width, calculated by the two 
methods , is approximately 1.5®. 

The antenna beam width, as defined by equation (4-1) , has been 
used for sensitivity calculation and will be used in the numerical 
calculations discussed in section 5. 

4.2 SIDE LOBES 

From figures 4-1, 4-2, and 4-3, it appears that the contribution 
of the power returned from the side lobes is not significant 
compared to that returned from the main lobe for 0 > 20®, but 
might be significant for 0 < 20®. To carefully examine the 
relative levels of the main and side lobes, the two-way gains at 
incidence angles 0® and 10® are plotted as a function of the 
crosstrack angle <() in figures 4-6 and 4-7 for HH and W polariza- 
tions. For W polarization, the first side lobe levels are 
approximately 11 dB and 14 dB down from the main lobe at 0 = 0® 
and 0 = 10®, respectively. On the other hand, the HH polarization 
first side lobe levels are approximately 27 dB and 32 dB down from 
the main lobe at 0 = 0® and 0 = 10®, respectively. It is also 
clear that the size of the first side lobe for HH polarization is 
much smaller than for W polarization. Therefore, it is expected 
that there are more side lobe contributions from W polarization 
than from HH polarization. The gain of each of the remaining 
side lobes is comparable to that of the first side lobe in each 
of the two polarizations. In general, the side lobe gains become 
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Figure 4-6.— The variation of the two-way antenna 
gains with crosstrack angles. 
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progressively smaller as |<}»1 increases. For the estimation of 
the side lobe effect, it is desirable to obtain the total contri- 
bution from all of the side lobes. 


The derivation of the backscattering coefficient a® normally 
includes the ground cell associated only with the main lobe an- 
tenna pattern. However, the power return at the scatterometer 
receiver includes contributions from other areas, where the side 
lobe gains are finite, and from the antenna coupling effect, 
which is discussed in section 4.3. At a given Doppler frequency 
band, associated with the alongtrack incidence angle 0 and 
angular width A0, the side lobe contribution to the power return is 


PgO) 


/ 




G.G a® 
t r 


dA 


(4-2) 


where the integration is carried out over the area defined by the 
Doppler frequency band and all side lobes. To make a rough but 
reasonable numerical estimate of Pg(0), the area element of inte- 
gration is to be approximated by the one shown in figure 4-8. 

The two isodopplers determine the total ground area, which back 

scatters the RP power in the corresponding Doppler frequency band. 

• 

Curves and L 2 represent the surface ]oci of the 3-dB points of 
a given side lobe. Point is the intersection between the iso- 
doppler, corresponding to the central Doppler frequency and the 
curve tracing the two-way gain maximum of the side lobe. The area 
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AA = Ax • Ay 


Figure 4-8.— An area element used in the calculation of 
side lobe contribution to total signal. 
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element from this side lobe is approximately given by: 


= AX^AY^ 


(4-3) 


and equation (4-2) can be written: 


Ps(6) 




N 


(4ir) 


E 

i=l 


G.a 


1 1 
T~ 


AX.AY. 


(4-4) 


2 

where is the side lobe maximum two-way gain at point for 
the area element AA^. The values of a? and refer to point 
Cj^. The summation is over all the side lobe area elements within 
the two isodopplers. 


The estimated values for P_(0) are presented in section 5, where 

s 

the results from all numerical calculations are summarized and 
discussed. 


4,3 ANTENNA CROSS-POLARIZATION AND ITS COUPLING EFFECT 
Since the nature of the horizontal and vertical antennas and the 
isolation between them are not perfect, undesired signals are 
generated because of antenna coupling effects which appear in the 
receiver data channel as portions of the backscattered signals. 

The relative contribution of the undesired signals to the total 
signal power level has important bearing in the data analysis and 
in the interpretation of the backscattering process of the RF waves. 
Thus, it is necessary to estimate the relative signal contribution 
resulting from the effect of antenna coupling in ordev to correctly 
interpret the measured results. 
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Since the backscattered signal in the 0.4-GHz scatterometer is 
delay gated with a predetermined time period for corresponding 
aircraft altitude which ensures that the return signals are 
only from the target, the direct power leakage into the receiver 
data channels results from the coupling between the transmit/ 
receive antennas would be gated out and would not appear as power 
return in the data channel. The main source of the undesired 
signals are from the coupling effect of the antenna cross-polari- 
zation of the horizontal and vertical antennas. Antenna cross- 
polarizations are mainly due to the nature of the practical 
antenna, where a horizontal antenna will generate, simultaneously, 
a vertically polarized pattern in a much smaller magnitude when 
compared to the main horizontally polarized pattern; the same 
is true for a vertical antenna. 


To see how this source of undesired signals are generated, a 
specific example is described below. Note that the expression 
PrHH^®)' which is the power level of the signal return at an 
incidence angle 0 in the horizontal-transmit and horizontal- 
receive data channel including the antenna coupling effect, can 
be written by; 


X^AAP 




( 470 ^( 6 ) 


(6) L 


G‘{e)a«H(9) + 


+ Gj^(6)GHc(e)aj^(e) + G^^(0)a^{e) 


] 


(4-5) 
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where 


Gjj = one-way horizontal antenna gain, either in the transmitting 
or the receiving mode of operktibn 
Grc = one-way cross-polarized horizontal antenna gain, either in 
the transmitting or the receiving mode of operation 

The first term inside the bracket in equation (4-5) is the domi- 
nant term which represents the backscattered signals of interest. 
The three remaining terms are contributions of the undesired 
signals resulting from antenna cross-polarization and, thus, 
its coupling effect. These three terms in sequential order, 
represent: 

a. The return signal in the horizontal receiver caused by the 
signal transmitted by the cross-polarized component of the 
horizontal antenna, backscattered in the horizontal polari- 
zation, and received by the prime component of the horizon- 
tal antenna 

b. The return signal in the horizontal receiver caused by the 
signal transmitted by the prime component of the horizontal 
antenna, backscattered in the vertical polarization, and 
received by the cross-polarized component of the horizontal 
antenna 

c. The return signal in the horizontal receiver caused by the 
signal transmitted by the cross-polarized component of the 
horizontal antenna, backscattered in the vertical polarization, 
and received by the cross-polarized component of the horizontal 
antenna 
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Expressions similar to equation (4-5) were derived for the other 
three power returns of different polarizations, namely, ^rVH' 

and The results were summarized in table 4-2, including 

the one for discussed above. For each of the four power 

returns in the table, the first term always gives the dominant 
term of interest from which the associated backscattering coeffi- 
cient a® is derived. The other three terms, representing the 
undesired signal returns, would have to be estimated and their 
relative contribution to the total power level would have to be 
assessed in the derivation of a®. In table 4-2, if the cross- 
polarized backscattering coefficients o®^^^ and cf®^^ are assumed 
to be the same (namely, incidence angle) , 

then further simplification of the results are given in table 
4-3. The assumption is reasonable because backscattering coeffi- 
cients obey the reciprocity theorem in general, and the experimen- 
tally measured data of the two coefficients usually show that 
they are about the same in magnitude. Table 4-3 will be used in 
the calculations presented in section 5. 

To show the level of antenna coupling, the peak one-way gains of 
both prime and cross (coupling) polarizations for the horizontal 
and vertical antennas were plotted as a function of 0 in the aft 
direction from 0° to 70® (figs. 4-9 and 4-10) . For the horizon- 
tal antenna (fig. 4-9) , the cross-polarization level is generally 
about 16 dB or more below that of the prime polarization; it 
reaches a maximum difference of 25 dB at 0 = 22®. To reduce the 
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Figure 4~9.— The variation of peak one-way gain of the horizontal 
antenna with incidence angle for both prime and cross 
polarizations. 
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Figure 4-10.— The variation of peak one-way gain of the vertical 
antenna with incidence angle for both prime and cross 
polarizations . 
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TABLE 4-2.- THE COMPONENTS OF THT RACKSCATTERED SIGNAL 
IN RECEIVER DATA CHANNELS 


Power return 

Components of 
backscattered signal 

P 

1. 

^H^HH^H 

*^rHH 




2. 



3. 

G CT® ,G 

. - - 


H HV'^HC 


4. 

^HC^^W^HC 

P 

1. 

^h'^hvS 

*^rHV 



2. 

^HC^W^V 


3. 

Vhh'Vc 


4. 

^HC^^S^C 


1. 


P 



rVH 




2. 



3. 



4, 




■ysiiymsiHHH 

P 

1. 


rW 




2. 

*^VC^HV^V 

■ 

3. 


■' 

4. 

^vc^hhSc 
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TABLE 4-3.- THE COMPONENTS OF THE BACKSCATTERED SIGNAL 
IN RECEIVER DATA CHANNELS 


[Assuming 
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coupling effect, the backseat ter ing coefficient should be measured 
in the region from 15® to 35® where the cross polarization level 
is about 20 dB below the prime polarization level. 

For the vertical antenna (fig. 4-10) , the cross-polarization level 
is generally about 11 dB or more below that of the prime polariza- 
tion level and attains a maximum difference of 13 dB at 
0 = 40® to 44®. Comparing figures 4-9 and 4-10, it is clear that 
the horizontal antenna has less antenna coupling effect than the 
vertical antenna, especially for 0 in the range of 15® to 35®. 

Generally, the peak gain cross-polarization level (figs. 4-9 and 
4-10) is not located at <(> = 0®, where the peak prime polarization 
gain usually would be located. This is illustrated in figures 
4-11, 4-12, 4-13, and 4-14, where the one-way gains of both prime 
and cross polarizations were plotted as a function of at 0 = 30® 
and 0 = 60® in the aft direction for the horizontal and vertical 
antennas. The peak cross-polarization gain occurs at ({) = -10® 
for both the horizontal and vertical antennas when the incidence 
angle is 30®, and it occurs at ip = -25® when the incidence angle 
becomes 60®. Also, from figures 4-11 through 4-14, it is seen 
that at 0 = 60® the cross-polarization level is higher than at 
0 = 30® for both the horizontal and vertical antennas. This is 
in agreement with what is shown in figures 4-9 and 4-10. There 
are two cross-polarization lobes in the antenna main beam (figs. 
4-12, 4-13, and 4-14); only one with comparatively smaller 
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Fiqure 4-11.— The variation of one-way horizontal antenna gain 
with crosstrack angle for both prime and cross polarizations 
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Figure 4-13.— The variation of one-way vertical antenna gain with 
crosstrack angle for both prime and cross polarizations. 
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Figure 4-14.- The variation of one-way vertical antenna gain with 
crosstrack angle for both prime and cross polarizations. 


4-28 



magnitude appears in figure 4-11. This indicates that the 
horizontal antenna at 0 = 30® contains less antenna coupling 
effects when compared to the other three cases shown in figures 
4-12 through 4-14. 

Based on figures 4-11, 4-12, 4-13, and 4-14 and table 4-3, the 
relative contributions caused by antenna coupling could be esti- 
mated, at some selected locations, for each of the four power 
returns P^^^, and P^^. 

The results are summarized in table 4-4 for 0 = 30® and 0 = 60°, 
each with two values of n> where the prime and cross-polarization 
gains attain the peak values. For 0 = 30®, the values for (j) are 
0® and -10®; for 0 = 60®, 0® and -25°, respectively. All of the 
relative power returns and a°'s in the table are expressed in dB. 

The components of the backscattered signals are listed in the 
same sequence as those in table 4-3. Since the values of 0 °^^^^ 
and 0 °^ are generally observed to be about 10 to 15 dB greater 
than those for or °®^^ (ref. 14), it is immediately clear 

that for sod the contributions from the undesired 

signals are approximately 20 dB or more below that of the signal 
of interest at 0 = 30®. At 0 = 60° and (|> = -25°, the difference 
between the signal of interest and the undesired signals becomes 
less than 10 dB for but is about 20 dB for For the 


4-29 



TABLE 4-4.- THE RELATIVE CONTRIBUTIONS TO THE 
TOTAL BACKSCATTERED SIGNALS 



Components of backscattered signal, 

dB 1 

Power return 

0 = 30“ 

0 = 60“ 1 


(J) = 0“ 

^ -10“ 

c{) = 0“ 

<P = -25“ I 

^^rHH 


+ 8.4 

a“nH + 28.8 



o“hv - 13.1 

-10.6 

"“hv - 



°°HV - 

-10.6 

“°HV- 



a<>vv - 38.0 

-29.6 

- ^2.0 


^^rHV 

^**HV ^ 15.7 

+10.6 

o“hv + 24.2 

+13.4 

(or AP VH) 
r 

- ^-2 

- 8.4 

a“vv “ 6.2 

+ 4.4 


^‘’hh ■ 

+ 1.2 

^“hh 

+ 8.8 . 


q“hv “30.5 

-17.8 

°°HV - 24-2 

- 0.2 


'’"vv 

+12.8 

°°w * 

+ 7.6 


"°HV - 

+ 3.4 




“'hV - 

+ 3.4 




a=HH - 23.0 

- 6.0 

“"hh - I®-" 


Antenna gains 

Gjj = 5.9 

4.2 

G„ = 14.4 
11 

9.6 


Ghc = -19.0 

-14.8 

Chc = -16.0 

0.6 


Gv = 9-8 

6.4 

G^ = 9.8 

3.8 


Sc = “^^'5 

-3.0 

Gyc = -8.2 

-0.8 
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remaining two power returns and the undesired signal 

contributions may be significant or even larger than the signal 
of interest. Yor example, at <|> = -10“ and 0 = 30®, the third 
term undesired signal level may be as high as the signal of 
interest; whereas at (j> = -25® and 0 = 60®, it may reach approxi- 
mately 5 dB larger than the signal of interest. These results 
suggest that, depending upon locations, the contributions of the 
undesired signal could be three times higher than the signal of 
interest for either ^Pj-jjy ^^rVH* 

The power at a given Doppler frequency band, as observed by the 
0.4 -GHz scatterometer , represents the total backscattered signal 
from the entire ground cell corresponding to that frequency band. 
Therefore, the comparisons made in table 4-4 at selected locations 
are clearly not sufficient to give a valuable assessment of the 
antenna coupling effect. A reasonable approach for estimating 
the effect is to integrate the power of each of the backscattered 
signal components, desired and undesired, over the entire ground 
cell and compare the results. This is done by numerical inte- 
gration described in section 5. 

It is noted (fig. 4-11) that the measurements of the antenna gains 
are limited to the levels above approximately -28 dB. This 
suggests that the actual cross-polarized gain levels could be 
lower than the measured values in some regions. However, the 
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estimate on the total undesired signal contribution is heavily 
weighted towards the components having the higher power returns 
associated with higher antenna cross-polarization levels. Con- 
sequently, the saturation in the cross-polarized gain measurements 
does not introduce significant errors in the numerical estimates 
of the antenna coupling effect performed in section 5. 

The coupling in the side lobes also makes undesired contributions 
to the total signal power return. However, the contribution to 
the total signal caused by the prime side lobe gains is small 
(sections 4.2 and 5.3). The undesired signal level caused by the 
antenna coupling in the side lobes is expected to be even smaller; 
therefore, these should not be explored. 
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5. NUMERICAL RESULTS AND DISCUSSION 


5.1 INTRODUCTION 

The numerical computations of the receiver signal output for 
various polarizations were performed according to equation (3-22) . 
The area elements with the relative coordinates are shown in 
figure 5-1. The center of the area element lies on the isodoppler 
corresponding to the incidence angle 0^ and the crosstrack angle 
4)^. The integration of is made from alongtrack (4)^^ = 0) to 
offtrack ( 4 )^^ ^ 0) depending on the antenna beam width which must 
be included. The offtrack incidence angle 0^ is related to the 
alongtrack (y-axis) incidence angle 0^ by: 

0. = cos ^ (cos 0 cos 4> . ) (5-1) 

1 o 1 


and the offtrack distance to alongtrack distance R(®q) hy: 

R(0^) = R(0^)/cos 4>i 

(5-2) 

= h/cos 0^ 

The area element AA^ is given by: 

AAj^ = Ax^ . Ay^ (5-3) 

where Ay^ for offtrack AA^ is related to alongtrack Ay^ by: 
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Figure 5-1.— The coordinates of the area element AA^. 


(5-4) 


AYi = Ay^cos 4)^ 

where Ay^ is 50 meters (164 feet) in the computation. The width 
Ax^ is defined to be the distance substended by 2° of the cross- 
track angular increment and can be written as: 


Ax . 
1 


h 


tan((J)^ + 1) - tan (4)^ 



(5-5) 


where h is the aircraft altitude and is taken to be 460 meters 
(1510 feet) . 


The antenna gains G^(0^) and G^(G^), the backscattering coeffi- 
cient a®(0^),and the distance R(0^) inside the summation sign 
of equation (3-22) were evaluated at the centers of the AA^ ' s . 

When these centers did not fall on the positions where the actual 
measurements of ^^(0^) and G^(G^) were made, linearly interpolated 
values from the immediate neighbors were used. The summation was 
carried out to 4>^ = ±60°, where the antenna crosstrack patterns 
indicate that almost all the powers are inside this angular range. 
To make the estimate of the antenna coupling effect simple and 
reasonable, a single ±60° crosstrack boundary was applied to the 
calculations of all four components of desired and undesired 
signal returns for the three polarizations as shown in table 4-3. 

The computer program, designed and implemented by J. R. Wang 
(ref. 10) for the 1.6-GHz scatterometer , was modified to calculate 
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the receiver power return for the 0.4-GHz scatterometer . The 
modified version given in the appendix can evaluate the power out- 
puts in accordance with equation (3-22) . Eight angles of incidence 
at 5°, 10°, 15°, 20°, 30°, 40°, 50°, and 60° were chosen for the 
computations. Since the alongtrack cell length is a fixed value 
of 50 meters (164 feet) , the bandwidth Af^^ changes when the angle 
of incidence decreases from 22 Hz at 0 = 5° to 3.8 Hz at 0 = 60°. 

The values of the antenna feed-line loss L and the calibration 

s 

constant K were treated as constants in the program. 

The program, coded in Fortran, consisted of one main program and 
two subroutines. The first subroutine was used to obtain the 
values of the two-way antenna gain for the calculation of signal 
returns, given the incidence angle 0^ and the crosstrack angle 
(|)^. The second subroutine was used to evaluate the backscattering 
coefficient o°, given the incidence angle 0^ at the area element 
AA^ . The functional dependence of and 0 °^^^ on 0^ is given 

in equations (3-27) and (3-28) for water and land, respectively. 

The values of a°„„ (or a°,,„) were assumed to be 15 dB lower than 
a°uu or a°^_, at all incidence angles. The program was executed 

nn VV 

on the UNIVAC 1110 computer at JSC. 

The two-way antenna gains for each of 10 combinations in table 
4-3 were previously saved on a magnetic tape, one file for each 
combination. For a a° versus 0 curve, all the signal outputs 
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corresponding to 10 components listed in table 4-3 can be obtained 
with a single computer run. 

5.2 THE 0 ° DEVIATION 

As discussed in section 3, it is necessary to assume a functional 
dependence of o® on the incidence angle 0 in equation (3-22) of 
a® derivation. This assumption was made according to equations 
(3-27) and (3-28) for backscattering from water and land, respec- 
tively. The receiver signal outputs in terms of (B/A)/Afp from 
the computer runs are shown in figures 5-2 and 5-3 for backscat- 
tering over water and land, respectively. These signal outputs 
could be compared with those actually obtained from observations 
to determine the validity of the assumed a® dependence on G. 

The procedure is straightforward and will not be performed here. 
Instead, a comparison will be made between the first [as repre- 
sented by equation (3-21) ] and the second [as represented by 
equation (3-22)] methods of derivation. To make this comparison, 
the data displayed in figures 5-2 and 5-3 were taken to be the 
receiver output signals for both methods. Then the a® dependence 
on 0 , as derived from the second method, will be given by equation 
(3-27) for water and equation (3-28) for land. The numerical 
values from these equations at the incidence angles of 5®, 10®, 
15®, 20°, 30°, 40°, 50®, and 60® are listed under the o° columns 
in tables 5-1 and 5-2. The o°'s for both linear and cross polari- 
zations over water and land were given in tables 5-1 and 5-2, 
respectively. Note that no distinction is made between o°„„ 

nn 
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COMPUTED BACKSCATTERED POWER (B/A)/ f^^AdB 


10 


+ 

X 





+ 

X 


X 


o 
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X 

o 
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o 

X 




O HH POLARIZATION 

XHV (OR VH) POLARIZATION 

+ VV POLARIZATION 

1 I I I I I 

0 10 20 30 40 50 60 70 

INCIDENCE ANGLE 6. DEGREE 

Figure 5-3.- The computed backscattered signals over 
land for all four polarizations. 
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and and between and in this second method of a° 

derivation (section 5.1). 

The corresponding values of the o^'s derived from the first 
method of equation (3-21) were given in tables 5-1 and 5-2 under 
the ^°HV same numerical values for 

the parameters such as altitude, aircraft speed, cell length, 
antenna feed loss, and calibration constant K were used in 
equations (3-21) and (3-22). From tables 5-1 and 5-2, it is seen 
that the results of the linear polarizations are more similar for 
the two method than are the results of cross polarization power 
return. This is because the cross polarization power return 
contains significant amounts of the antenna coupling term, which 
will be discussed in later subsections. The inclusion of the 
antenna coupling effect increases the results derived from the 
first method by 2 to 5 dB higher than from the second method. 

For linear polarization, with negligible antenna coupling effects 
the results derived from the first method is approximately 0 to 
2 dB lower than from the second method. 

5.3 EFFECT OF ANTENNA SIDE LOBES 

The side lobe contribution to the total power return can be esti- 
mated by using equation (4-4) . However, to compare the power 
levels between the main beam and the side lobe signal returns, 
it is more convenient to write equation (4-4) in a form similar 
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TABLE 5-1. - A COMPARISON OF o“ VALUES 
DERIVED FROM TWO DIFFERENT APPROACHES OVER WATER 


[dB] 


Incidence 

ann 1 p 0 . 

Linear polarization 

Cross polarization 

degree 

a° 


^“vv 


^‘’hv 

5 

- 7.57 

- 9.80 

- 9.68 

-22.57 

-20.36 

10 

-12.00 

-13.30 

-13.49 

-27.00 

-24.01 

15 

-16.14 

-17.37 

-17.83 

-31.14 

-28.38 

20 

-20.00 

-22.20 

-21.74 

-35.00 

-33.14 

30 

-26.85 

-28.80 

-28.17 

-41.85 

-40.48 

40 

-32.56 

-33.95 

-33.10 

-47.56 

-45.45 

50 

-37.14 

-37.85 

-36.69 

-52.14 

-48.64 

60 

-40.56 

-40.74 

-39.30 

-55.56 

-50.07 
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TABLE 5-2. - A COMPARISON OF o® VALUES DERIVED 
FROM TWO DIFFERENT APPROACHES OVER LAND 


[dB] 


Incidence 
angle 0/ 
degree 

Linear polarization 

Cross polarization 

a® 

^°HH 

°°W 

0° 


5 

5.01 

3.02 

3.14 

- 9.99 

- 7.53 

10 

0.59 

- 0.45 

- 0.63 

-14.41 

-11.12 

15 

- 3.40 

- 4.35 

- 4.82 

-18.40 

-15.33 

20 

- 6.97 

- 8.90 

- 8.45 

-21.97 

-19.81 

30 

-12.85 

-14.56 

-13.94 

-27.85 

-26.20 

40 

-17.05 

-18.21 

-17.38 

-32.05 

-29.66 

50 

-19.55 

-20.08 

-18.93 

-34.55 

-30.76 

60 

-20.37 

-20.33 

-18,94 

-35.37 

-29.49 
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to equation (3-22) . This is readily done and the final expression 
for the side lobe signal return is: 


(B/A)/ fj^ 


(47t) ^ K Af, 


N 


E 



Ax. 

1 



(5-6) 


where N is the number of the side lobe area elements contributing 

to the power return at the Doppler frequency f^ and bandwidth Af^; 

L and K have the same meaning and values as in equation (3-22) ; 
s 

2 

and G./ o° , Ax., and Ay. were defined in section 4.2. The varia- 
tion of a° with incidence angle G were again determined by equa- 
tions (3-27) and (3-28) for backscattering over water and land, 
respectively. From equation (5-6) , the side lobe power level can 
be estimated in terms of {B/A) /Af^ and compared directly with the 
main beam power return shown in figures 5-2 and 5-3. 


Only the W polarization contains a significant amount of side 
lobe power; therefore, only the contributions of power return 
from the first side lObe on both sides of the main beam of the 
W polarization were computed and plotted in figures 5-4 and 5-5 
for water and land, respectively. Also included in figures 5-4 
and 5-5 is the depressed 30 dB main lobe power return. From 
figures 5-4 and 5-5, it is clear that the side lobe power level 
for W polarization is about 30 dB or more lower than that of 
the main lobe power. As observed in section 4, the other polari- 
zations contain much smaller amounts of side lobe level as com- 
pared to W polarization; thus, it can be concluded that the side 
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ESTIMATED SIDE LOBE BACKSCATTERED POWER (B/Aj/Afp,, dB 



0 10 20 30 


INCIDENCE ANGLE 6, DEGREE 

Figure 5-4.— Estimated side lobe signal return over water 
corresponding to the indicated main beam incidence 
angle. 



ESTIMATED SIDE LOBE BACKSCATTERED POWER (B/A)/Af., dB 


-20 


+ SIDE LOBE SIGNAL (VV POLARIZATION) 
© MAIN LOBE SIGNAL (DEPRESSED 30 dB) 


INCIDENCE ANGLE 6. DEGREE 

Figure 5-5.— Estimated side lobe signal return over land 
corresponding to the indicated main beam incidence 
angle. 
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lobe contributions to the backscattered signal powers are insig- 
nificant for all polarizations. 


5 . 4 ANTENNA COUPLING EFFECT 

The computer outputs of the backscattered signals are displayed 
in sequential order in tables 5-3, 5-4, and 5-5 for HH, HV, and 
VV polarizations, respectively. The computed signal outputs were 
derived for backscattering over water and land and for incidence 
angles of 5°, 10°, 20°, 30°, 40°, 50°, and 60°. The assignment 
of the backscattered signal components under columns 1, 2, 3, 
and 4 was the same as that given in table 4-3. Consequently, 
the values in column 1 were the backscattered power of interest, 
and those under the other three columns were the backscattered 
powers caused by antenna coupling. The values in the Total 
column give the sum of the first four columns. The values for 
G°,„, or o°,,„ were taken to be 15 dB less than those for o° 

nV Vn nil 

and in the computations. It is clear from tables 5-3 and 

5-5 that by comparing columns 1 and Total, the contribution to 
the backscattered signal caused by antenna coupling is totally 
nejligible in these cases. Therefore, it is not necessary to 
consider the effect of antenna coupling when deriving or ^°yy 

The results shown in table 5-4 indicate a significant contribution 
from columns 2 and 3 to the total power. The contributions from 
column 2 aie larger than those of column 1, which indicate that 
the antenna coupling effect overshadows the power return of interest. 

P^GFi IS 


original 
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TABLE 5-3.- NUMERICAL VALUES OF THE BACKSCATTERED SIGNALS FOR 
HORIZONTAL TRANSMIT AND HORIZONTAL RECEIVE 


Incidence 

angle G, _ 

degree 1 


15-33 

18.58 

21.40 

22.85 

20.33 

17.09 

15.76 

16.24 


HH (wat( 
3 


•49.20 

•52.55 

•56.78 

•60.99 

-58.66 

•52.72 

-49.55 

-48.21 


-49. 

, 20 

-52. 

,55 

-56. 

, 78 

-60. 

,99 

-58. 

, 66 

-52. 

. 72 

-49. 

. 55 

-48. 

, 21 


■64.25 
•53 .85 
■48.18 
■44.49 


[IN TERMS OF B/A/Af^] 


Total 


-2.510 
-5.723 
-8.384 
-9 . 549 
-6.088 
-1.357 
2.018 
4.162 


•18.57 
■21.40 
■22.85 
■20 .33 
•17.09 
5 

■16.23 


-36.36 
-39.64 
-43.68 
-47.63 
-44 .36 
-36.84 
-31.52 


■36 .36 
■39.64 
■43.68 
■47.63 
•44 .36 
■36.84 
■31.52 
■27.41 


■35.19 
■38.97 
■45.18 
■53.48 
■49.82 
■37 . 74 
•29.78 
•23.23 
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TABLE 5-4.- NUMERICAL VALUES OF THE BACKSCATTERED SIGNALS FOR 
HORIZONTAL TRANSMIT AND VERTICAL RECEIVE 


[IN TERMS OF B/A/Af^] 


Incidence 
angle 0 , 
degree 


HV or 

VH (water) 


HV or VH (land) 

1 

2 

3 

4 

Total 

1 

2 

3 

4 

Total 

5 

-14.05 

-12.63 

-17.59 

— 

-43.22 

- 9.53 

-1.237 

0.2031 

-4 . 742 

-30.36 

3.297 

10 

-17.21 

-15.51 

-21.13 

-46.22 

-12.61 

-4.358 

-2.617 

-8.210 

-33.27 

0.275 

15 

-19.71 

-17.90 

-25.42 

-50.39 

-15.26 

-6.698 

-4.837 

-12 .33 

-37.23 

-2.213 

20 

-20.61 

-19.34 

-29.08 

-55.60 

-16.66 

-7.324 

-5.985 

-15.72 

-42.13 

-3.334 

30 

-18.51 

-18.70 

-27.55 

-54.26 

-15.32 

-4.277 

-4.361 

-13.28 

-39.83 

-1.041 

40 

-17.18 

-17.58 

-23.87 

-49.50 

-13.90 

-1.462 

-1.723 

-8.042 

-33.45 

1.887 

50 

-18.55 

-17.99 

-23.23 

-47.87 

-14.61 

-0 . 786 

-0.0275 

-5.266 

-29.58 

3.27 

60 

-22.45 

-19.75 

-25.32 

-47.70 

-17.16 

-2.079 

0.8982 

-4.618 

-26.55 

3.418 


TABLE 5-5.- NUMERICAL VALUES OF THE BACKSCATTERED SIGNALS FOR 
VERTICAL TRANSMIT AND VERTICAL RECEIVE 


[IN TERMS OF B/A/Af^^] 


VV (water) 


W (land) 


■ 8.414 

■11.47 

-13.63 

■13.92 

•12.25 

■12.81 

•16.88 

■24.23 


2 

3 

4 

Total 

1 

2 

— 

3 

— 

4 

Total 

-36.76 

-36.76 

-33.03 

-8 . 386 

4 . 406 

-23.92 

-23.92 

-20.16 

4.434 

-39.71 

-39.71 

-35.54 

-11.44 

1.392 

-26 . 81 

-26 . 81 

-22.58 

1.422 

-42.20 

-42.20 

-36.99 

-13.60 

-0.619 

-29.13 

-29 . 13 

-23.85 

-0.586 

-43.50 

-43.50 

-38.05 

-13.89 

-0.631 

-30.15 

-30.15 

-24.61 

-0.604 

-43.27 

-43.27 

-39.23 

-12.23 

1.986 

-28.95 

-28.95 

-24.79 

2.002 

-43.95 

-43.95 

-40.36 

-12.80 

2.907 

-28.12 

-28.12 

-24 . 36 

2.922 

-46.92 

-46.92 

-42.64 

-16.86 

0.870 

-29.01 

-29.01 

-24.45 

0.892 

-52.06 

-52.06 

-45.60 

-24.18 

-3.870 

-31.48 

-31.48 

-24.62 

-3.819 





While there were comparable contributions at all eight incidence 

angles under investigation, the least severe ones appeared to be 

at 0 = 30® where the value of column 2 is about the same as that 

of column 1. At all other incidence angles, the values of column 

2 are larger than those of column 1. To account for the antenna 

coupling effect in the derivation of a®„,, or o°,,rr, a correction 

factor shall be included in the data processing algorithm. The 

values of the correction factor will be the difference between the 

columns labelled Total and 1 in table 5-4, assuming or a®,,„ 

to be 15 dB less than those of and 

nn VV 

The effect of the antenna coupling on the backscattered signal 
power level is expected to depend strongly on the relative values 
of or with respect to and cr°yy To make a rough 

estimate of this dependence, a similar computation was performed 
with and again given by equations (3-27) and (3-28) for 

backscatter ing over water and land, respectively. But the values 
of o®j^y or were assumed to be 10 dB less than those of o®j^pj 

and The results of this computation are shown in tables 

5-6, 5-7, and 5-8 for HH, HV, and W polarizations, respectively. 
Again, it is clear from tables 5-6 and 5-8 that the contribution 
to the backscattered signal caused by antenna coupling is still 
negligible in the linear polarization cases. Therefore, it is not 
necessary to consider the effect of antenna coupling when a°„,, was 
assumed to be either 15 dB or 10 dB less than those for o®u„ for 

HH 

the linear polarization cases. 
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TABLE 5-6.- NUMERICAL VALUES OF THE BACKSCATTERED SIGNALS FOR 
HORIZONTAL TRANSMIT AND HORIZONTAL RECEIVE 


[IN TERMS OF B/A/Af^,] 


Incidence 
angle 0 , 
degree 

HH (water) 

— . . _ . 

HH (land) 

1 

2 

3 

4 

Total 

1 

2 

3 

4 

Total 

5 

-15.33 

-44.20 

-44.20 

-48.08 

-15.32 

-2.510 

-31.36 

-31.36 

-35.19 

-2.496 

10 

-18.58 

-47.55 

-47.55 

-51.95 

-18.57 

-5.723 

-34.64 

-34.64 

-38.97 

-5.710 

15 

-21.40 

-51.78 

-51.78 

-58.37 

-21.39 

-8.384 

-36.88 

-38.68 

-45.18 

-8.375 

20 

-22.85 

-55.99 

-55.99 

-66.96 

-22.85 

-9.549 

-42.63 

-42.63 

-53.48 

-9.545 

30 

-20.33 

-53.66 

-53.66 

-64.25 

-20.33 

-6.088 

-39.36 

-39.36 

-49.82 

-6.084 

40 

-17.09 

-47.72 

-47.72 

-53.85 

-17.08 

-1.357 

-31.84 

-31.84 

-37.74 

-1.348 

50 

-15.76 

-44.55 

-44.55 

-48.18 

-15.75 

-2.018 

-26.52 

-26.52 

-29.78 

-1.980 

60 

-16.24 

-43.21 

-43.21 

-44.49 

-16.22 

4.162 

-22.41 

-22.41 

-23.23 

4.189 












TABLE 5-7.- NUMERICAL VALUES OF THE BACKSCATTERED SIGNALS FOR 
HORIZONTAL TRANSMIT AND VERTICAL RECEIVE 


[IN TERMS OF B/A/Af^] 


Incidence 
angle G, 
degree 


HV 

or VH (water) 


HV or VH (land) 

1 

2 

3 

4 

Total 

1 

2 

3 

4 

Total 

5 

- 9.05 

-12.63 

-17.59 

-38.22 

- 7.06 

3.763 

0.2031 

-4.742 

-25.36 

5 . 757 

10 

-12.21 

-15.51 

-21.13 

-41.22 

-10.18 

0.642 

-2.617 

-8.210 

-28.27 

2.693 

15 

-14.71 

-17.90 

-25.42 

-45.39 

-12.76 

-1.698 

-4.837 

-12.33 

-32.23 

0.269 

20 

-15.61 

-19.34 

-29.08 

-50.60 

-13.94 

-2.324 

-5.985 

-15.72 

-37.13 

-0.632 

30 

-13.51 

-18.70 

-27.55 

-49.26 

-12.23 

0.723 

-4.361 

-13.28 

-34.83 

2.027 

40 

-12.18 

-17.58 

-23.87 

-44.50 

-10.86 

3.538 

-1.723 

-8.042 

-28.45 

4.899 

50 

-13.55 

-17.99 

-23.23 

-42.87 

-11.88 

4.214 

-0.0275 

-5.266 

-24.58 

5.948 

60 

-17.45 

-19.75 

-25.32 

-42 . 70 

-15.01 

2.921 

0,-8982 

-4.618 

-21.55 

. ___ ■ 

5.492 






TABLE 5-8.- NUMERICAL VALUES OF THE BACKSCATTERED SIGNALS FOR 
VERTICAL TRANSMIT AND VERTICAL RECEIVE 


[IN TERMS OF B/A/Af^^] 


Incidence 
angle 0, 
degree 

W (water) 

... 

W (land) 

1 

2 

3 

4 

Total 

1 

2 

3 

4 

Total 

5 

-8.414 

-31.76 

-31.76 

-33.03 

-8.359 

4.406 

-18.92 

-18.92 

-20.16 

4.461 

10 

-11.47 

-34.71 

-34.71 

-35.54 

-11.41 

1.392 

-21.81 

-21.81 

-22.58 

1.451 

15 

-13.63 

-37.20 

-37.20 

-36.99 

-13.57 

-0.619 

-24.13 

-24.13 

-23.85 

-0.560 

20 

-13.92 

-38.50 

-38.50 

-38.05 

-13.87 

-0.631 

-25.15 

-25.15 

-24.61 

-0.583 

30 

-12.25 

-38.27 

-38.27 

-39.23 

-12.22 

1.986 

-23.95 

-23.95 

-24.79 

2.017 

40 

-12.81 

-38.95 

-38.95 

-40.36 

-12.78 

2.907 

-23.12 

-23.12 

-24.36 

2.937 

50 

-16.88 

-41.92 

-41.92 

-42.64 

-16.84 

0.870 

-24.01 

-24.01 

-24.45 

0.911 

60 

-24.23 

-47.06 

-47.06 

-45^60 

-24.15 

-3.870 

-26.48 

-26.48 

-24.62 

-3.787 


It can be shown that even if the values of or o°,,„ are com- 

rlV VM 

parable to those of and the effect of antenna coupling 

remains negligible in and cr®^ estimations. For both water 

and land backscatterings , the results in table 5-7 indicate a 
1- to 2.5-dB contribution caused by antenna coupling, which is 
smaller than that shown in table 5-4. The contributions from 
column 2 become smaller than those of column 1 because of increas- 
ing values of o° . . 


The smallest coupling effect occurs at^ 0 = 30®, where the values 
in column 2 are about 5 dB less than those in column 1. By com- 
paring tables 5-4 and 5-7, it can be seen that the larger the 
difference is between the assumed values of a®„,, and a°„„, the 
larger the contributions are from the antenna coupling effect. 
This is because a decrease in o®^^^ results in a decrease in 
column 1; simultaneously, the unchanged values in column 2 will 
result in a larger share of contribution from column 2. This is 
what happens in column 1 of tables 5-4 and 5-7, where exactly 
5 dB difference is shown between the two tables. As a result 
of the undesirable contribution from antenna coupling, or 

a°vH would be overestimated in the data processing algorithm. 

A correction would be required which would depend on the relative 
level of a®„,, with respect to o®„„. For o®,.,, to be 10 dB less 

than those of a°„„ and a®..,,, the value difference between the 

rln VV 

Total column and colum.n 1 in table 5-7 should be employed. 


5- 22 



5.5 OUTPUT SIGNAL PRECISION 


5.5.1 ANTENNA PATTERN MEASUREMENTS 

The uncertainty in the antenna pattern measurements was not 
recorded and supplied with the gain values. Since the New Mexico 
State University measured both the 1.6-GHz and 0.4-GHz antenna 
patterns, it is assumed that the same degree of uncertainty is 
placed to the 0.4-GHz system. An uncertainty of approximately 
±0.5 to ± 1 dB in one-way gain results in the uncertainty of 
± 0.7 to ± 1.3 dB for the two-way gain used in the a° estimates. 
Another source of uncertainty, which is significant for the 0.4- 
GHz scat ceroivieter because of its large antenna array, comes from 
the fact that the gain measurements were made without a mockup 
to simulate the effect of the C-130 aircraft. It is expected 
that, after the antennas were installed in the aircraft, the 
gain pattern would differ from that measured in the laboratory. 
The magnitude of this antenna pattern change cannot be estimated. 

5.5.2 MEASUREMENT OF CALIBRATION CONSTANTS 

An uncertainty exists in the measurements of the calibration 
constants K for various polarizations and gain settings. A 
reasonable estimate placed this uncertainty to be no more than 
±0.5 dB for each of these measurements. 

5.5.3 AIRCRAFT PARAMETERS 

Factors associated with aircraft motion, such as roll, pitch, 
drift, vertical velocity, speed, and altitude changes, could 
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introduce errors in the received signals. The effects of roll, 
pitch, and vertical velocity alter the incidence angle 0 and were 
studied by Krishen (ref. 4) with respect to the 13.3-GHz scatter- 
ometer system. Following the same approach, the errors introduced 
by these factors in the o° estimates are presented in tables 
5-9, 5-10, and 5-11 for all four polarizations. In making these 
error calculations, the aircraft speed is assumed to be 77 m/sec 
(253 ft/sec), the vertical velocity is 4 m/sec (13.1 ft/sec), and 
both roll and pitch angles are 5°. Clearly, the change in 0 , 
caused by pitch and vertical speed, could be either plus or minus. 
Tables 5-10 and 5-11 list larger modified incidence angles only. 

Changes in both aircraft speed and drift would introduce changes 
in the Doppler frequency which is associated with a given G. For 
example, at 0 = 30°, a change of speed from 77 m/sec (253 ft/sec) 
to 70 m/sec (230 ft/sec) would give a corresponding change in 
Doppler frequency from 103 Hz to 93 Hz. Similarly, a drift of 
approximately 10° would shift the Doppler frequency from 103 Hz 
to 101 Hz. These Doppler frequency changes, if not corrected for 
in data processing, could result in erroneous o° estimates. 

The variation in aircraft altitude would also introduce error in 
the a° estimate. From equation (3-22), o° is observed to be pro- 
portional to and inversely proportional to the area A of the 

2 

ground cell. Since R is proportional to altitude h and A to h , 

2 

0 ° varies with h . Thus, if h changes from 460 to 440 meters 
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TABLE 5-9. - ERRORS INI -ODUCED IN THE o° 
ESTIMATE DUE TO AN AIRCRAFT ROLL OF 5® 


Incidence 
angle, i 

degree 

Modified 

incidence 

angle, 

degree 

Error in a°(dB) 

HH 

HV or VH 


5 

7.07 

-0.25 

-0.35 


10 

11.17 

-0.29 

-0.41 


15 

15.79 

-0.40 

-0.40 


20 

20.59 

-0.44 

-0.53 


30 

30.38 

-0.38 

-0.38 


40 

40.26 

-0.13 

-0.09 


50 

50.18 

I 

-0.05 

0 



W 


-0.50 

-0.59 

-0.40 

-0.59 

-0.38 

-0.07 

0.05 


60 


60.13 


- 0.02 


0.05 


0.07 









TABLE 5-10. - ERRORS INTRODUCED IN THE o° ESTIMATE 


DUE TO AN AIRCRAFT PITCH OF 5® 


Incidence 

angle, 

degree 

Modified 

incidence 

angle, 

degree 

Error in o° (dB) 

HH 

HV or VH 


5 

10 

-0.25 

-0.75 

-1.25 

10 

15 

-1.50 

-2.00 

-2.50 

15 

20 

-3.50 

1 

to 

o 

-4.00 

20 

25 

-4.50 

-4.65 

-4.75 

30 

35 

-5.00 

-3.95 

-3.25 

40 

45 

-3.00 

-1.75 

-0.50 

50 

55 

-1.25 

0.15 

1.50 

60 

65 

-0.50 



1.65 

3.75 
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TABLE 5-11. - ERRORS INTRODUCED IN THE o® ESTIMATE 
DUE TO AN AIRCRAFT VERTICAL SPEED OF 4 M/SEC 

(13 ft/sec) 



Incidence 

angle, 

degree 

Modified 

incidence 

Error in a° (dB) 

angle , 
degree 

HH 

HV or VH 

W 

5 

7.97 

-0.25 

-0.45 

-0.75 

10 

12.97 

-0.75 

-1.10 

-1.50 

15 

17.97 

-1.50 

-2.00 

-2.50 

20 

22.97 

-2.50 

-2.80 

-3.00 

30 

32.97 

-3.00 

-2.60 

-2.50 

40 

42.97 

-1.75 

-1.10 

-0.50 

50 

52.97 

-0.75 

0 

0.75 

60 

62.97 

-0.50 

0.80 

2.00 












(1510 to 1443 feet) , the error in the o° estimate would be 
approximately 0.39 dB. 

5.5.4 COMMENT 

The uncertainty introduced in the antenna pattern measurements 
and by the presence of the aircraft body is difficult to correct. 
But this uncertainty can be treated as a bias and should not 
affect the repeatability of the data taken from one mission at 
another level. The calibrate signal level could be adjusted to 
accommodate a particular mission. For each adjustment, a new 
calibration constant K would have to be measured; therefore, 
the measurement should be made as precisel\ as possible. The 
changes in aircraft parameters could introduce large error in 
0 ° estimation (tables 5-9, 5-10, 5-11). This type of error could 
be minimized by making corrections to the aircraft parameters 
during data processing. The error sources from the antenna 
coupling could be corrected by proper data reduction and manipu- 
lation using tables 5-4 and 5-7 for cress polarization power 
returns. However, if the measured data indicate an irregular 
variation, the data should be used with reservation. 
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6 . 


CONCLUSIONS AND RECOMMENDATIONS 


A system analysis was performed on the 0.4-GHz airborne scattero- 
meter. The principle of the scatterometer operation, calibration 
procedures, resolution cell length requirement, and data handling 
were described. The antenna gain patterns were studit 1, and the 
effect of the antenna cross coupling effects were analyzed. 

Several error sources in the estimation of a® were briefly dis- 
cussed. In the course of the entire analysis, only the immediate 
relationship between the power input to the power output from 
the scatterometer system was examined. The interacting process 
of the RF waves with terrain surfaces was completely omitted. 

The terrain surface characteristics were reflected through the 
assumed functional forms of the backscattering coefficient a® (G) . 
The major conclusions that resulted from this study are listed 
below: 

a. To obtain a sufficient independent sample within a resolution 
cell, the cell length should be made larger than 75 m as des- 
cribed in section 3.3, especially at large incidence angles. 

b. The gain patterns for the vertical and horizontal antenna are 
quite different. The gain maximum occurs at 0 = 66° for the 
horizontal antenna, and it occurs at 0 = 42® for the vertical 
antenna. The beam width variations (with 0) of the two 
antennas closely follow one another, with the vertical antenna 
having the smaller beam width. 
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c. The cross-polarization effects, measured for the two antennas, 

have important bearing on the receiver signal outputs of both 
HV and VH polarizations. If or o°.,„ is less than 

rlV Vn htn 

and hy approximately 15 dB, the power contributions 

caused by the antenna cross-polarization effects were larger 
than the return power of interest for the HV or VH receiver. 

If or is less than and <^®yv approximately 

10 dB, the cross polarization effects were found to be com- 
parable to the return power of interest. Thus, in general 
conditions, a correction factor should be added to make the 
estimation meaningful for the HV or VH receiver. The magni- 
tude of the correction may be obtained from tables 5-4 and 
5-7. 

d. The side lobe contributions to the receiver signal outputs 
are negligible. For the cases examined, the backscattered 
power levels from the side lobes are approximately 30 dB or 
more below that from the main beam. 

e. The variations of the aircraft parameters affect the estima- 
tion of a“. For example, a pitch of 5° could result in 

an error of approximately 5.0 dB in some cases. 

The following recommendations should be consjderud jn order to 
improve the performance of the system and the quality of the 
final report. 

a. The effect of the aircraft body on the antenna pattern cannot 
■ be estimated. To improve the accur<^c\ of o° measurements. 
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it is desirable to calibrate the aircraft scatterometer with 


a well-calibrated ground system on a homogeneous terrain 
surface . 

b. The antenna cross-polarization effects should be considered 

in the estimations of or If information provided 

by o° or o° is indispensable in the aircraft remote 
sensing of soil moisture, carefully interpolated correction 
factors should be used to obtain a reasonable estimate of 

or a° . In other cases, the measured data of or 

nV Vn nV 

should be used with reservation. 

Vn 

c. since aircraft parameter variations significantly affect the 
estimation of a°, corrections to the variation of the air- 
craft parameters should be made in the data processing. 
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APPENDIX 


THE COMPUTER PROGRAM LISTING 


1 ♦ 


DIMENSION ThETA(O),P0*(8t2),BUF(92J 

2 * 


DATA TMtTA/5.,10.,15.,2G.,30.,40,,5Q., 

3 * 


COMMON/li/IDADO,I)(,IY,lXMAX,IYMAX 

4 ♦ 

C 

EME‘^ CONSTANT PARAMETERS. 

5* 


AL0SS:8 .4 

6 ♦ 


PI tr 3. 14 1 6 

7* 


HA =4 6 U. 

8 * 


VA = 77 . 

9 * 


SL AM U = 0 .75 

in* 


CELEN^Sa. 

11 • 


PI :P IE/ 180. 

12* 


IF IL -1 

13* 


TX =92 

14* 


lY =36 

15 * 


I X P A X ; 9 

16 * 


I Y ma X : 3b 

17* 


CALL PIN I T ( I DADD , NURO S > 

18 * 


IE0F=IFIL*10 

19* 


50 CONTINUl 

20 * 

C 

GET DISC AREA 

21 * 


\S =1 F IL 

22 * 


K = ID A DO 

23 * 

c 

PUT ANTENNA FOOTOPINT ON DISC 

24 * 


lOG 9EAD(l,40Q0,END = 2CO)(dUF(NR),NP = l ,92) 

25 * 

c 

WRITE(6,5000)(3UF(NR),NP=1,92) 

26* 


call RwR ITE ( H,3UF , IX , NSTAT ) 

27* 

5 

IFtNSTAT.EO.l) GO TO 5 

28* 


K = K I X 

29* 


GO TO luO 

30* 


200 CONTINUE 

31 * 

c 

ENTER *-AIN LOOP OF COMPUTATION. 

32* 

c 

alonu-tpack COhduTaTICN. 

33 * 


00 530 1=1,8 

34* 


TMTA=Th£TA( I )*PI 

35* 


FD0P-2.*VA*SIN(THTA)/SLAMD 

36* 


DF=2.*VA*(C0S(THTA)**3)*CELEN/HA/SLAMD 

37* 


CL AMD = SL« mO**2 

38 * 


CP 1= <4.*PIE ) **3 

39* 


SUM= [: . 

40 * 

c 

cross-track COmDJTATION. 

4 1 * 


I F ( I . L T . 7 ) J J= 2 1 

42 * 


IF (I .GF .7) JJ=21 

4 3 * 


no 6 no M = 1, 2 

44 * 


no 4nr j=i,jj 

45 * 


f’Hl: ( J-l ) *PI*2. 

46 * 


PH IL - ( ( J- 1 ) *2. - i . ) *PI 

47 * 


PH TP = ( ( j- 1 ) *2. ♦ 1 . ) *PI 

48 * 


OX = H A * f T A N( P HIP ) -T A N ( PH I L ) ) 

49 * 


ny =cELrN/cos (ph i ) 

50 * 


OA =DX *DY 

51 * 


d=HA/COS<THTA)/COSIPHI) 

5? * 


PP zR * *4 , 

S3* 


TH INC=AC0S( ha/p ) 

54 * 


TH INC^THINC /pi 

55* 


CAIL t>COFF( T HINC , SI GMA , N S, M ) 

56 * 


AI=ATAN«TAN(ThTA)/COS(PHT) ) 

57* 


AT =A I/PI 

58* 


CALL GAT (PGL ,PGR ,I ,J, AI) 

59* 


SI G= 10.** (SI DMA /ID.) 

60* 


IF (J . EO. 1 )PGP = G , 

61 • 


SUM;Sum*( (PGL*P bP)*SIG*OA/PR) 


60 . / 


A-1 


62 * 

C 

wRITE(6,bC0D)P5L,PGR,SIGMA,SIb,nA,RK,SUM 

63 * 

4on 

CONT INUE 

64 * 

C CALCULATE POWER OUTPUT EUR AN INCIDENT ANGLE. 

65* 


PA C- SUM# C la MO/C PI /DF 

6 6 * 


PAC= in.#AL0G10( PAC ) 

67* 


IF(NS.Eg.I.OP.NS.E0.7.OR.NS.E0.d) SCAl=122.3 

68 * • ' 


IF(NS.EU.2.0R.NS.E0.3.0R.NS.E0.4.0R.NS.EO.‘=) 

69 * 


IF (NS. FO, 6. OR. NS. EC. 9. OR. NS. EQ. 10) SCAL-123. 

70 * 


POi*( I,M)rPAC-AL05S*SCAL 

71 * 

6C0 

CONT INUE 

72 * 

C 

WRITE(6,6COO)PAC,EnOP,DE,CLAHD,CPI,ALOSS,SCAL 

73 * 

5 DC 

CONT INUE 

74 « 

C 

WR IT E (6 , 3'IOD ) IE IL , NS , M 



wR IT E (6 , 1 COD ) 

7b * 


wRITE(6,2Q0D)((P0W(I,M),I=l,8),Mrl.2) 

77* 


IF IL - IFIL-^I 

78* 


IFdFIL.LT. IE OF) GO TO 5C 

79* 


STOP 

80 * 

lonu 

F0PMAT(2H , ’INCIDENT ANGLE A T * , 3 X , • 5 ’ , 1 1 X , ’ 1 C 

81 * 

c 

13x,’30’,10X’4Q»iGX,’50’,inx,’60*,/) 

82* 

2D0C 

60PM A T ( 16 X, 8 D12 .** ) 

83 * 

3GCD 

FORM A T ( 3 I 10 ) 

84 * 

4 COO 

FORMAT (22E6 . 1 ) 

85* 

50UC 

F0RMAT(20F6. 1) 

86 * 

6QnC 

F0PMAT(7D12.4) 

87* 


FN P 


SCAL: 13S.7 


' ,10X , M5* 


IJX 


•23’ 


1 <■ 
? <« 

3 * 

4 * 

5 » 

6 ■<' 

7 • 

8 * 
9* 

n * 

1 I * 
12 * 
13 * 
14 ♦ 
15* 
16* 
17* 
18 * 
IQ# 

23 • 

21 * 


SJ t-R CuT INf 3AT ( GL ,GR , I , J ,A I > 

CO^'MON./ci/lDAnD,IX,IY,lXMflX,IYMAX 
dimension BJFK92) ,BUF2(92) 

c determine The PvOSITion of The antemna footprint. 
fli -a 1/2. 

TA::3 6- AINT(AT) 

AA 1= 36. -A I 
I API - I A- 1 
IP 1 = I A* iX ♦! DADO 
IP2- ( IA-1 )* IX4-TDADD 
C READ antenna PECURDS. 

call RPLAD( IPl.BUF 1, IX,I ST ) 

5 IF(I5T.EC.1)00 TO 5 
IF(IST,NF.3)WRITF(6, IDOOIA 
CALL KRlAD( IR2.BDF2, IX, I ST ) 

6 tf(ist.lc.1)Go to 6 

IF(IST.NE.D)WRITF(6, ILGOIAPI 
C PERFORM linear INTERPOLATION. 

IF(j.Ea.l)GC TO 7 
PLUPUFI(47-J) 

PL 2-:6oF2 ( 47-JJ 


A- 2 


2 ?* 

23 * 

24 * 
2S* 
26 « 
27* 
28 * 
29 * 
3H* 

31 * 

32 * 

33 * 

34 * 
35* 
36* 
37* 

38 * 

39 * 
4D * 
4 1* 

42 * 

43 * 

44 * 
4 5* 
46 * 
47* 
4 8* 


1 * 

? * 

3 * 

4 * 

5 * 

6 * 


7 


8 

C 

laoG 

2 onc 

3onc 

4noo 


P9 1 = I^UF 1 ( 45* J) 
PQ2- 5UF2 ( 45* J) 


PL IL : 10. **( P LI/ n. ) 

TF lJ.6T.5.A\D.PLl.5w.b. )PL 1L=Q, 

PL2L-10.**( »L2/1L'. > 
IF(J.6T.5.ANn.PL?.EJ.o.)PL2L=0. 

PR 1L=1D.**( PRl/in. ) 
IF(J.GT,5.AND.PR1 .LQ.c.)9R1L:0. 

PR2L : 10. ♦*( PP2/ 1 n. ) 
TF(J.GT.5.AND.PR2.FC.C.)PR2L;0. 

GL = ( ( (AAI-IA )/( lAPl-lA) )*(Pl 2L-PL1L ) )*PLlL 
GRH((AAI-IA)/(IAP1-IA))*(PR2L-PR1 L))^pR1L 
50 TO 8 
PI =BUF1 (46) 

P2 :BUF2 (4b) 

P1-1G.**{P1/10. ) 

P2=1C.**(P2/ 10. ) 

GL:( ( (AAI-IA )/( lAPl-IA) )*(P2*PD) *P1 

GR :Gl 


CONTI NUt 

WRITL(6,4u0D)P1»P2,PL1,PL2,PR1,PR2,GL,6R,AAI 

RETURN 

FORMAKIH ,’DISC READ PROBLEM FOR A\6(.E*,I5) 
FORM AT ( 1 H , F 6.2 ,5 X , 3 I 10 ) 

FORMAT (IH , 7F12 .4 ) 

POPMAT (IH , 9F8.2 ,2110) 

EN D 


I A, TAPI 


su f-.R nuT I N L 
IF (M . EO. 1 > 
I F ( M . f 0 . 2 ) 
I F ( N . L 0 . 1 . 
RETURN 
END 


eCQFF (A ,SIC*-A,N,M) 

SIGHArC.0"bn*A**Sfl,971*A-2.86 
EIGMArG.0G844*A**?-l,ni*A^9.85 
OR.N.EC.2.0R.N.E0.5.0R.N.L0.6) SI5 


‘'ArSTGHA-l 


OHV 

< )V' 
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